{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# ``atmosphere.py`` tutorial\n",
    "\n",
    "This tutorial needs your help to make it better!\n",
    "\n",
    "This tutorial requires pvlib > 0.6.0.\n",
    "\n",
    "Authors:\n",
    "* Will Holmgren (@wholmgren), University of Arizona. 2015, March 2016, August 2018."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "# built in python modules\n",
    "import datetime\n",
    "import logging\n",
    "import os\n",
    "import inspect\n",
    "\n",
    "# python add-ons\n",
    "import numpy as np\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pvlib\n",
    "from pvlib.location import Location"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "tus = Location(32.2, -111, 'US/Arizona', 700, 'Tucson')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Location: \n",
      "  name: Tucson\n",
      "  latitude: 32.2\n",
      "  longitude: -111\n",
      "  altitude: 700\n",
      "  tz: US/Arizona\n"
     ]
    }
   ],
   "source": [
    "print(tus)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                           apparent_zenith      zenith  apparent_elevation  \\\n",
      "2014-06-24 00:00:00-07:00       124.039086  124.039086          -34.039086   \n",
      "2014-06-24 00:01:00-07:00       124.065488  124.065488          -34.065488   \n",
      "2014-06-24 00:02:00-07:00       124.090882  124.090882          -34.090882   \n",
      "2014-06-24 00:03:00-07:00       124.115267  124.115267          -34.115267   \n",
      "2014-06-24 00:04:00-07:00       124.138641  124.138641          -34.138641   \n",
      "\n",
      "                           elevation     azimuth  equation_of_time  \n",
      "2014-06-24 00:00:00-07:00 -34.039086  352.695507         -2.366772  \n",
      "2014-06-24 00:01:00-07:00 -34.065488  352.970528         -2.366922  \n",
      "2014-06-24 00:02:00-07:00 -34.090882  353.245771         -2.367072  \n",
      "2014-06-24 00:03:00-07:00 -34.115267  353.521228         -2.367222  \n",
      "2014-06-24 00:04:00-07:00 -34.138641  353.796891         -2.367372  \n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOzdeVzVVf748de57Pu+g4LKIggiIqKWoeaSpWZTqdWMTU3WVE7NzLeyfu1fK2scZ75m1lgpNZppmkXmlLmFO4rgAiKgrIKyL7Jz7/n9ATIuiKDgZTnPx4MH8Lmf5X2v+L7nns857yOklCiKoii9k0bfASiKoihdRyV5RVGUXkwleUVRlF5MJXlFUZReTCV5RVGUXkwleUVRlF7MUN8BXMnR0VF6e3vrOwxFUZQeJT4+vkhK6XTl9m6X5L29vTl8+LC+w1AURelRhBBZrW1X3TWKoii9mEryiqIovZhK8oqiKL1Yt+uTb01DQwO5ubnU1tbqOxRFz0xNTfH09MTIyEjfoShKj9Ajknxubi5WVlZ4e3sjhNB3OIqeSCkpLi4mNzcXHx8ffYejKD1Cj+iuqa2txcHBQSX4Pk4IgYODg/pEpygd0COSPNDhBC91kqryOgpzKqmraeyiqJRbTb3RK0rH9Ijumo6QUlJf08iF0jq0jTqAlu+Koih9TYda8kIIUyFEnBDiqBAiSQjxVvP2N4UQZ4UQic1fUy855mUhRLoQ4pQQYnJnP4FLNTZoKS+oobywBgRYO5p15eW6nXffffeWXCcmJoZFixYB8N1335GcnNzyWFRUlJrMpig3oKsWcOpod00dMF5KORQIBaYIISKbH/uHlDK0+WsLgBAiEJgNBAFTgOVCCINOiv0yVeV1lORX0VCnxdLOBHs3C4xNmy/Vwxa/amy8se6lW5Xkp0+fzoIFC4Crk7yiKB2XdaKYdQvjKDtf3enn7lB3jWx6q7nQ/KtR81dbKXQG8LWUsg7IEEKkAxHA/huIFYC3fkgiOa/iqu06rUTqJAaGGrik27ahTouBoQaNwbX7cgPdrXljWtB1r33vvfeSk5NDbW0tzz33HPPmzcPS0pInn3ySnTt3Ymdnx9dff42TkxNRUVGEhoYSFxdHRUUFK1euJCIigri4OJ5//nlqamowMzNj1apV+Pv7Ex0dzY8//khtbS1VVVX88MMPzJ8/n+PHj9PY2Mibb77JjBkziI6OJiYmhurqak6fPs3MmTP54IMPWLBgATU1NYSGhhIUFMSaNWuuiv+TTz7hk08+AaC8vBxvb2927tzJ1q1beeONN6irq2PgwIGsWrUKS0tLvL29mTt3Lj/88AMNDQ188803BAQEEB0dzeHDh3nooYeIiYnh119/ZeHChWzcuBGAb775hqeffpqysjI+//xzbr/99uu+torSF5Wdr2bPhjSyjhdj62JObXVDp1+jwzdehRAGQohEoAD4RUp5sPmhZ4UQx4QQK4UQds3bPICcSw7Pbd525TnnCSEOCyEOFxYWdjQkADQGAgOjyxN8Z1u5ciXx8fEcPnyYpUuXUlxcTFVVFWFhYRw5coQ77riDt956q2X/qqoq9u3bx/Lly3nssccACAgIIDY2loSEBN5++21eeeWVlv3379/PF198wY4dO3jnnXcYP348hw4dYufOnbzwwgtUVVUBkJiYyLp16zh+/Djr1q0jJyeHRYsWYWZmRmJiYqsJHuCpp54iMTGRQ4cO4enpyV/+8heKiopYuHAh27Zt48iRI4SHh7NkyZKWYxwdHTly5Ah//OMfWbx48WXnGz16NNOnT+dvf/sbiYmJDBw4EGj6JBIXF8c///nPy14PRVGa1FU3sHdDGmvfPkheWhmjfzOI2a9F4Opj0+nX6vCNVymlFggVQtgCm4QQQ4CPgf+lqVX/v8DfgcdoPeVe1fKXUq4AVgCEh4e32bnSnhb3RTqtjqLcC1jamWJubdzu465l6dKlbNq0CYCcnBzS0tLQaDTMmjULgEceeYT77ruvZf85c+YAMHbsWCoqKigrK6OyspK5c+eSlpaGEIKGhv++c0+cOBF7e3sAtm7dSkxMTEtira2tJTs7G4AJEyZgY9P0xxAYGEhWVhZeXl7tfh7PPfcc48ePZ9q0aWzevJnk5GTGjBkDQH19PaNGjWrZ9+LzGT58ON9++227zn/pMZmZme2OS1F6O61WR1JsHoc2Z1Bb3cDgUW6MnDEACxuTLrvmDY+ukVKWCSF2AVOklC1NPCHEp8Dm5l9zgUuzjyeQd6PX1Kddu3axbds29u/fj7m5OVFRUa2O1750iN+Vw/2EELz22muMGzeOTZs2kZmZSVRUVMvjFhYWLT9LKdm4cSP+/v6XnePgwYOYmPz3D8LAwKBDffjR0dFkZWWxbNmylutMnDiRtWvXtrr/xWt15Do3coyi9GZSSrKOF7N3Yzpl56vx8LdjzP2DcPKy6vJrd3R0jVNzCx4hhBlwJ5AihHC7ZLeZwInmn2OA2UIIEyGED+ALxN182LdeeXk5dnZ2mJubk5KSwoEDBwDQ6XRs2LABgK+++orbbrut5Zh169YBsGfPHmxsbLCxsaG8vBwPj6Yeq+jo6Gteb/LkyXz44Yctd9wTEhKuG6ORkdFlnwyuFB8fz+LFi1m9ejUaTdM/fWRkJHv37iU9PR2A6upqUlNTr3uti6ysrKisrGz3/orS1xTlVhLzf4n8uPwYAFOfDmHG86G3JMFDx1vybsAXzSNkNMB6KeVmIcS/hRChNHXFZAJPAkgpk4QQ64FkoBF4prm7p8eZMmUKn3zyCSEhIfj7+xMZ2TSoyMLCgqSkJIYPH46NjU1LYgews7Nj9OjRLTdeAV588UXmzp3LkiVLGD9+/DWv99prr/H8888TEhKClBJvb282b958zf0B5s2bR0hICGFhYa32yy9btoySkhLGjRsHQHh4OJ999hnR0dHMmTOHuro6ABYuXIifn1+7XpfZs2fzxBNPsHTp0pY3O0VRoKK4hkObM0g5cA4Tc0Nun+VL0FgPDAxu7RxU0VVjM29UeHi4vHKc9cmTJxk8eHCHz9XZffKtsbS05MKFC1dtj4qKYvHixYSHh3fJdfuyG/17UJRbobqinsP/ySQp9ixCCIKjPBh+lzemFl1bVE8IES+lvCrh9LoZr4qiKPpQV91AwtZsju7IQdsoGTzGjRFTvbG0M9VrXCrJ36TWWvHQdKNWX4qLi5kwYcJV27dv346Dg4MeIlKU3qu+ppHjv+aSsDWbuupGfMOdiZg2AFsXc32HBqgk3ys5ODiQmJio7zAUpVerrWrg2M5cju3Ioa66kf5DHBg5Y8Atu6HaXirJK4qidEB1RT1Ht+dw/NdcGmq1+Ax1JHyqN879rfUdWqtUklcURWmHiqIaju3IJWn3WRobdQwa7szwKd44elrqO7Q2qSSvKIrShnNnyknclsOZhAIQAr8IF4ZP6Y+dq8X1D+4GVJK/RWJiYkhOTmbBggV89913+Pn5ERgYCKjhlorS3ei0Os4kFpG4LZvzGRUYmxkSOrEfwVGeWNnrd7RMR/WRJK//uQDTp09n+vTpQFN53nvuuaclySuK0j1UV9Rzcl8eSbF5VJbUYu1kxu2zfAkY5Yaxac9Mlz1m+T99+uSTTwgNDSU0NBQfHx/GjRvH1q1bGTVqFGFhYTzwwAMtQym9vb154403CAsLIzg4mJSUFKCphMGzzz7Lvn37iImJ4YUXXiA0NJTTp08DTeV5IyIi8PPzY/fu3Xp7rorS10gpyU0p4acVJ/hiwV4OfHcGaydT7noqmIffiiRknFePTfDQE1vy/1kA5463a1eBxLa5njxtTSV2DYa7Fl3z4aeeeoqnnnqKhoYGxo8fz2OPPdZSntfCwoL333+fJUuW8PrrrwP/Lc+7fPlyFi9ezGeffdZyrovlee+55x7uv//+lu0Xy/Nu2bKFt956i23btrXrOSqKcmOqK+o5dfAcSbvPUl5Qg4mFIcHjPQm6zb3H9Le3R89L8np0sUSvnZ2dKs+rKD1QY4OWjKNFnDp4juykEqRO4jbIhhF3+zAwzAlDoy5ZuE6vel6Sb6PFfSWpk5TlVGJpZ4K59c3Va760RO+PP/6oyvMqSg8hpST/dDmnDpwjPb6A+ppGLGxNGDaxH/4jXbF37z2t9tb0vCSvBxdL9O7evRuNRkNkZCTPPPMM6enpDBo0iOrqanJzc9tduVGV51WUriWlpCCrktPxBaQfKaCyuBZDYw0DhznjH+mKh78dGk0XLiPXjXQoyQshTIFYwKT52A1SyjeEEPbAOsCbplLDD0opS5uPeRl4HNACf5JS/txp0d8irZXoVeV5FaV7kVJSkFlJ+pECTjcndo1G4DnYjoh7fBgwzKlH30C9UR0qNSyaljqykFJeEEIYAXuA54D7gBIp5SIhxALATkr5khAiEFhL0+Ld7sA2wK+tmvKdWmpYJynqpO4apftQpYaVi3RaHfmny8k4VsSZI4VUltSiMRB4DbZnYJgzPkMdu7zEb3fRKaWGZdM7wsWyi0bNXxKYAUQ1b/8C2AW81Lz9ayllHZAhhEinKeHv7/hTUBRFaSoMlp1cTOaxYrKTiqmrbkRj2JTYI6b54B3SdxJ7e3T4s0vzqlDxwCDgIynlQSGEi5QyH0BKmS+EcG7e3QM4cMnhuc3bFEVR2kVKSdn5ajKPF5N5rIj80+VIncTMygifoY54hzjiNdi+T3bFtEeHX5XmrpbQ5rVeNwkhhrSxe2t3Nq7qHxJCzAPmAfTr16+jIV2f/ie8KorSAbVVDeSmlJKTXEzOyVIqS2oBsHe3YNikfviEOOLsbd1nbp7ejBt+65NSlgkhdgFTgPNCCLfmVrwbUNC8Wy7gdclhnkBeK+daAayApj75G41JUZSeSduo43xGOdnJJeQkl1CQXQkSjE0N8PC3Y9ikfvQf4oC1o5m+Q+1xOjq6xgloaE7wZsCdwPtADDAXWNT8/fvmQ2KAr4QQS2i68eoLxHVS7Iqi9FAX118+m1pGXmopuallNNZpERqBi7c1I+72wWuwPS7eVmhu8cLXvU1HW/JuwBfN/fIaYL2UcrMQYj+wXgjxOJANPAAgpUwSQqwHkoFG4Jm2RtZ0FfXRQFH0S9ugoyCrgrz0MvJSy8g/XU5DXVMqsHEyI2CkK16B9nj422FipvrWO1NHR9ccA4a1sr0YuHpR0abH3gHeuaHoFEXpkRrqtZw/U05eWhl56WWcO1OBtkEHNPWr+490xd3PFvdBtljYquHNXUm9ZfYi7777Lq+88soNHevt7c3hw4dxdHTssnhGjx7Nvn37Ou38SvdRX9NI/unmpJ5WRkFWBTqtRAhw9LJiyO0euPvZ4jbIBjNLY32H26d0aDLUrdAVk6EsbE2wsOk5rYXGxkYMDTv+/mtpadlS8rijuiLJ30w8bVGTofSv9kJDU9dLc1IvyqlEStBoBM7eVrj72uI2qOlLdb/cGp0yGao7eD/ufVJKUtq9f0OtFo2haCo3fA0B9gG8FPHSdc917733kpOTQ21tLc899xzz5s3D0tKSJ598kp07d2JnZ8fXX3+Nk5MTUVFRhIaGEhcXR0VFBStXriQiIoK4uDief/55ampqMDMzY9WqVfj7+xMdHc2PP/5IbW0tVVVV/PDDD8yfP5/jx4/T2NjIm2++yYwZM4iOjiYmJobq6mpOnz7NzJkz+eCDD1iwYAE1NTWEhoYSFBTEmjVrWn0Oq1evZunSpdTX1zNy5EiWL1+OgYHBdfdZsWIFGRkZfPDBB0BTwbb4+Hg+/PDDVl+X1uK5mPSllLz44ov85z//QQjBq6++yqxZs9i1axdvvvkmjo6OnDhxguHDh7N69WqaJlor+lRVXkdeWhn5aWWcTSujJK8KAANDDS4+1gyf6o27ry2uPjYYmfS+So49WY9L8vq0cuVK7O3tqampYcSIEfzmN7+hqqqKsLAw/v73v/P222/z1ltvsWzZMgCqqqrYt28fsbGxPPbYY5w4cYKAgABiY2MxNDRk27ZtvPLKK2zcuBGA/fv3c+zYMezt7XnllVcYP348K1eupKysjIiICO68804AEhMTSUhIwMTEBH9/f+bPn8+iRYtYtmwZiYmJ14z/5MmTrFu3jr1792JkZMTTTz/NmjVr+N3vfnfdfe6//35GjRrVkuTXrVvH//t//++ar0tb8Xz77bckJiZy9OhRioqKGDFiBGPHjgUgISGBpKQk3N3dGTNmDHv37uW2227rhH89pSMqS2qbWumppeSll1N2vhoAQxMD3Aba4BvugrufLS79rTEwUqNfurMel+Tb0+K+qLO7a5YuXcqmTZsAyMnJIS0tDY1Gw6xZswB45JFHWurCA8yZMweAsWPHUlFRQVlZGZWVlcydO5e0tDSEEDQ0NLTsP3HiROzt7QHYunUrMTExLF68GIDa2lqys7MBmDBhAjY2NgAEBgaSlZWFl9el0xFat337duLj4xkxYgQANTU1ODs7t2sfJycnBgwYwIEDB/D19eXUqVMt9fRbe10cHByuGceePXuYM2cOBgYGuLi4cMcdd3Do0CGsra2JiIjA09MTgNDQUDIzM1WS72JSSiqKaslLK23pfqkoapp8ZGJuiNtAGwLHuOPua4tTP0s1pLGH6XFJviM680P+rl272LZtG/v378fc3JyoqChqa2uvvuYlXQtXdjMIIXjttdcYN24cmzZtIjMzk6ioqJbHLSz+W9daSsnGjRvx9/e/7BwHDx5sqT0PHas/L6Vk7ty5vPfeeze0z6xZs1i/fj0BAQHMnDkTIUS7X5crr3EtN/rclPaTUlJeWMPZU/9N6hdKm6qpmloY4e5rS8g4L9z9bHHwsFSzSns49ZbcTuXl5djZ2WFubk5KSgoHDjSV5NHpdC3lgr/66qvLWp3r1q0DmlquNjY22NjYUF5ejodHU/me6Ojoa15v8uTJfPjhhy0JMSEh4boxGhkZXfbJ4EoTJkxgw4YNFBQ0TUguKSkhKyur3fvcd999fPfdd6xdu7bl08u1Xpe24hk7dizr1q1Dq9VSWFhIbGwsERER131+yo2rrWogPb6AnatT+Per+1nz+gF2rTlFzskSXHxsGDvbj9mvRfDY327jrqeCGTrBCycvK5Xge4Fe3ZLvTFOmTOGTTz4hJCQEf39/IiMjgabWd1JSEsOHD8fGxqYlsQPY2dkxevTolhuvAC+++CJz585lyZIljB8//prXe+2113j++ecJCQlBSom3tzebN29uM8Z58+YREhJCWFhYqzdeAwMDWbhwIZMmTUKn02FkZMRHH31E//7927WPnZ0dgYGBJCcntyTla70ubcUzc+ZM9u/fz9ChQxFC8MEHH+Dq6tqy6Lly83Q6ybnT5WQnF7deJmBiPzwD7LB1MVc3tnu5Xj2EUuokhV08hPJawwSjoqJYvHgx4eFXjWhSbpIaQtm6+tpGck6WkHm0iMzjxdRWNbSUCfAKtFdlAnq5XjOEUlGU/6qvaeTM0ULSDxeQm1KKtlGHibkh/YMd8AlxwivQXo1T7+PUv/5NutZkn127dt3aQC5RXFzMhAlXV5nYvn17m6NelJ5B26AjK6mY1LjzZB4vQtugw8rBlCF3eOAT4ojrIBsMVGtdaaaSfC/k4ODQ5nh5pWcqL6wmKTaPk/vyqa1qwMzKiMAx7vhFuODiY6361pVWqSSvKN2YlJLck6UkbMsmJ7kEoRH4DHUk8DZ3vALsVP+6cl0qyStKNyR1koxjRcT/J5OCrEosbIyJmOZD4Bh3VbVR6ZCOLhriBXwJuAI6YIWU8v+EEG8CTwCFzbu+IqXc0nzMy8DjgBb4k5Ty506KvR0B37IrKUqnyU8vY883aRRkVWLtaErUw/4ERLqp8gHKDenoX00j8Fcp5WAgEnhGCBHY/Ng/pJShzV8XE3wgMBsIommZwOXNC470Ct7e3hQVFXXqOd99993Lfh89enSnnl/pvqrK6vj50xN8u/gIVWV1jP/dYB5+K5Kg2z1UglduWEcXDckH8pt/rhRCnAQ82jhkBvC1lLIOyBBCpAMRwP4bjLfXu7IGu6q/3vtJKUmNO8/udak0NugYcbc3wyb1V9UclU5xw80DIYQ3TatEHWze9KwQ4pgQYqUQwq55mweQc8lhubTypiCEmCeEOCyEOFxYWHjlw93C6tWriYiIIDQ0lCeffBKtVnvdxz/++GNefPHFln2io6OZP38+0FS2ePjw4QQFBbFixQqAy8rzPvzww0DTZCtoSgQvvPACQ4YMITg4uGVm7a5du4iKiuL+++8nICCAhx9+uM3aMEr3Ul/byNbPk9i2Khk7VwtmvxpBxLQBKsErneaGbrwKISyBjcDzUsoKIcTHwP/StJzq/wJ/Bx6j9V7xqzKQlHIFsAKaZry2de1z775L3cn2TX+XQEOdlguGos1xwyaDA3BtY0Wla5Xfvd7jqjyv0paygmr+88lxSvOriLx3AMMm9Ve1YpRO1+EkL4QwoinBr5FSfgsgpTx/yeOfAheLrOQCl9bA9QTybjhaPbleiV5VnlfpqMLsSn74MBGdTjJtfihegfb6DknppTo6ukYAnwMnpZRLLtnu1txfDzATONH8cwzwlRBiCeAO+AJxNxNwWy3uK0kpKcy++do11yq/e7GKpCrPq3TEuTPl/PDhUYxNDbjv+WHYupjrOySlF+ton/wY4LfAeCFEYvPXVOADIcRxIcQxYBzwZwApZRKwHkgGfgKekVJqr3Hubut6JXpVeV6lvYrzLrB52VHMrIy474XhKsErXa6jo2v20Ho/+5Y2jnkHeKeDcXWum7wPea3yu9d7XJXnVS5VVVbH5g+PYmCkYfqfQrGyN9V3SEof0LtLDV/srrExUbMEe5GeWGpYq9Xx/T8SKMy5wG9eCMPR00rfISm9zLVKDasZFopyCxz8/gz56eWMe9hfJXjlllJJXlG6WH56GQlbswm63R2/CFd9h6P0MSrJK0oX0jbo2Lk6BSt7U0b/ZpC+w1H6IJXkFaULJfySRem5au54yB9jU1X0Vbn1VJJXlC5SXVFP/M/ZDBzmRP8hakUuRT9UkleULnJ4SybaBh2R9w7UdyhKH6aSfBf4wx/+QHJycqec69LSw5mZmQwZMqRTzqt0rYriGpJ2nyVwjJua8KTolUryXeCzzz4jMDDw+ju2w5X15ZWe4ej2HJAw/C5vfYei9HG9Osl39sLGV5YHjomJITQ0lNDQUPz9/fHx8QEgKiqKixO6LC0teemllxg+fDh33nkncXFxREVFMWDAAGJiYoCmGjjPPvtsy3Xuuecedu3a1WrpYa1WyxNPPEFQUBCTJk2ipqamU5+jcvNqqxpI3pvPoBHOalaronc97nb/7vWpFOVcaPf+DXWNaAw0GBhe+/3M0cuS2x/0u+65riwP/Ouvv7aUBX7wwQe54447rjqmqqqKqKgo3n//fWbOnMmrr77KL7/8QnJyMnPnzmX69OnXvN6VpYczMzNJS0tj7dq1fPrppzz44INs3LiRRx555LqxK7dO0u6zNNZpGTaxn75DUZSel+T16VrlgT/44APMzMx45plnrjrG2NiYKVOmABAcHIyJiQlGRkYEBweTmZnZ4Rh8fHwIDQ0FYPjw4Td0DqXrSJ0kKTYPD387NbNV6RZ6XJJvT4v7UgVZFZjbmGB5k7VrrlUeePv27XzzzTfExsa2epyRkVFLt5FGo2kpC6zRaFpKAhsaGqLT6VqOaavs8JVlhVV3TfeSm1JKZUkto2aqETVK99ChPnkhhJcQYqcQ4qQQIkkI8VzzdnshxC9CiLTm73aXHPOyECJdCHFKCDG5s5/ArdJaeeCsrCyefvpp1q9fj5mZ2Q2f29vbm8TERHQ6HTk5OcTF/bfk/rVKDyvdU/LePEwsDPEJddR3KIoCdPzGayPwVynlYCASeEYIEQgsALZLKX2B7c2/0/zYbCAImAIsF0L0yMUrp0yZQmNjIyEhIbz22mtERkaSmZlJcXExM2fOJDQ0lKlTp97QuceMGYOPjw/BwcH8z//8D2FhYS2PXSw9fPHGq9J91Vyo50xiIf4jXTE06pF/5kovdFOlhoUQ3wPLmr+ipJT5Qgg3YJeU0l8I8TKAlPK95v1/Bt6UUu6/1jk7s9QwdF53jdJ9dNdSw8d35RL7dSqzXo3A0dNS3+EofUynlxoWQngDw4CDgMvF5f+av19cANUDyLnksNzmbYrS65xOKMTWxRwHDwt9h6IoLW4oyQshLGlazPt5KWVFW7u2su2qjw5CiHlCiMNCiMOFhYU3EpKi6FXNhXry0soYOMyp0+dnKMrN6HCSF0IY0ZTg10gpv23efL65m4bm7wXN23MBr0sO9wTyrjynlHKFlDJcShnu5OTU0ZAURe8yjhYhdZKBYc7X31lRbqGOjq4RwOfASSnlkkseigHmNv88F/j+ku2zhRAmQggfwBeIQ1F6mdNHCrB2NMXRS/XFK91LR8fJjwF+CxwXQiQ2b3sFWASsF0I8DmQDDwBIKZOEEOuBZJpG5jwjpdR2SuQd0r3WsVV6l7rqBnJTSgkZ76W6apRup0NJXkq5h9b72QEmXOOYd4B3OhiXovQYmceK0GklA4eprkal++nVBcoU5VY4nVCIha0JLt7W+g5FUa7SN5J8N++tKSsrY/ny5S2/5+Xlcf/999+Sa3/zzTcMHjyYcePGtfp4YmIiW7Zsafk9JiaGRYsW3ZLYeoL62kayk0uaRtVoVFeN0v30jSTfzV2Z5N3d3dmwYcMtufbnn3/O8uXL2blzZ6uPX5nkp0+fzoIFC25JbD1B1olitA06Boaprhqle+pxBcp2Rq+gIOtMu/dvqNWiMRAYGF37/cy5/wDGPTqvzfOsXr2apUuXUl9fz8iRI1m+fDlffvkl7733Hm5ubvj5+WFiYsKyZct49NFHueeee1pa45aWlly4cIELFy4wY8YMSktLaWhoYOHChcyYMYMFCxZw+vRpQkNDmThxIs888wz33HMPJ06coLa2lj/+8Y8cPnwYQ0NDlixZwrhx44iOjiYmJobq6mpOnz7NzJkz+eCDD64Z/9q1a3n33XeRUnL33Xfz/vvv8/bbb7Nnzx4yMjKYPn06f/vb3y47pr6+ntdff3IvHYUAACAASURBVJ2amhr27NnDyy+/TE1NDYcPH255nmZmZqSkpJCVlcWqVav44osv2L9/PyNHjiQ6OhqArVu38sYbb1BXV8fAgQNZtWoVlpa9YxTK6SOFmFkZ4TrQVt+hKEqrelyS14eTJ0+ybt069u7di5GREU8//TSrV6/mjTfeID4+HhsbG8aNG8ewYcPaPI+pqSmbNm3C2tqaoqIiIiMjmT59OosWLeLEiROX1Y2/6KOPPgLg+PHjpKSkMGnSJFJTU4GmVnZCQgImJib4+/szf/58vLy8rrpuXl4eL730EvHx8djZ2TFp0iS+++47Xn/9dXbs2MHixYsJD79qNjTGxsa8/fbbLUkdaEncF5WWlrJjxw5iYmKYNm0ae/fu5bPPPmPEiBEkJibi6enJwoUL2bZtGxYWFrz//vssWbKE119/vd2vf3fVWK8lK6kY/wgXNKqrRummelySv16L+0oFWRWYWxtjaXfjK/Rs376d+Ph4RowYAUBNTQ379u0jKiqKi5O3Zs2a1ZJ8r0VKySuvvEJsbCwajYazZ89y/vz5No/Zs2cP8+fPByAgIID+/fu3XGfChAnY2NgAEBgYSFZWVqtJ/tChQ5fF+vDDDxMbG8u9997bgVehddOmTUMIQXBwMC4uLgQHBwMQFBREZmYmubm5JCcnM2bMGKDp08GoUaNu+rrdQXZyCY11WgYOUxOglO6rxyV5fZBSMnfuXN57772Wbd99913LAiJXurQ+vJSS+vp6ANasWUNhYSHx8fEYGRnh7e3dZu34i8dfy5W15S/Wp+/IOW7WpfXxL43nYr18AwMDJk6cyNq1a7ssBn05nVCAiYUh7v6qq0bpvtSN13aYMGECGzZsoKCgqVpDSUkJw4YNY9euXRQXF9PQ0MA333zTsr+3tzfx8fEAfP/99y314MvLy3F2dsbIyIidO3eSlZUFgJWVFZWVla1ee+zYsaxZswaA1NRUsrOz8ff371D8I0eO5Ndff6WoqAitVsvatWtbXaqwNW3F1h6RkZHs3buX9PR0AKqrq6/7iacn0DboyDxWjM9QJwwM1H8jpfvq/X+dndBVGhgYyMKFC5k0aRIhISFMnDiR/Px83nzzTUaNGsWdd955WQ34J554gl9//ZWIiAgOHjyIhUVTVcKHH36Yw4cPEx4ezpo1awgICADAwcGBMWPGMGTIEF544YXLrv3000+j1WoJDg5m1qxZREdHX9Zibg83Nzfee+89xo0bx9ChQwkLC2PGjBntOnbcuHEkJycTGhrKunXrOnRdACcnJ6Kjo5kzZw4hISFERkaSkpLS4fN0NzknS6ivaVQToJRu76bqyXeFTq8nn12BudXN9cm3R3R09GU3KJWu0x3qyW+PTibjWBG//+C2NheJV5RbpdPryStKX6Vt1JFxrAifEEeV4JVurw/ceL01Q9seffRRHn300VtyrbaMHDmSurq6y7b9+9//bhn1ci0///wzL7300mXbfHx8rnlzuS/LOVlCXXUjA4erUTVK99cHknzfcvDgwRs6bvLkyUye3GPXWb+lTh8pwNjMEK8Ae32HoijX1WM+a3a3eweKfuj776CxXsuZxCJ8hjq2OYtaUbqLji4aslIIUSCEOHHJtjeFEGeFEInNX1MveexlIUS6EOKUEOKGm4mmpqYUFxfr/T+4ol9SSoqLizE17dqb6G05c7SQ+ppGAiJd9RaDonRER7trooFlwJdXbP+HlHLxpRuEEIHAbCAIcAe2CSH8bmTREE9PT3Jzc7mR9V8rS2oxNjXAxNyow8cq3Y+pqSmenp56u37K/nNY2Zvi4WentxgUpSM6umhIrBDCu527zwC+llLWARlCiHQgAtjfoQgBIyMjfHx8OnoYAJ/M30VIlCehvxl0Q8crykUXSmvJOVlC+FRvVVZY6TE6q1PxWSHEsebunItNHA8g55J9cpu3XUUIMU8IcVgIcfhGWuttUf8Vlc6Ssv8cSFRXjdKjdEaS/xgYCIQC+cDfm7e3ll9b7VSXUq6QUoZLKcMvFtFSlO5Eq9Vx4tdcPAPssHEy13c4itJuN53kpZTnpZRaKaUO+JSmLhloarlfWhLRE8i72evdCHW7VrlZp48UUFVez9AJV1f5VJTu7KaTvBDC7ZJfZwIXR97EALOFECZCCB/AF4i72et1mOqvUW6SlJKj23OxcTajf5CDvsNRlA7p0I1XIcRaIApwFELkAm8AUUKIUJoazJnAkwBSyiQhxHogGWgEnrmRkTWKom+5KaUUZFYwdrafuuGq9DgdHV0zp5XNn7ex/zvAOx0NSlG6CyklB2POYGlnwuAxbtc/QFG6md4/ZU8IUJOolBuUdaKY8xkVhE/1xtDIQN/hKEqH9f4kryg3SNuoY++GdKydzAgYrVrxSs+kkryiXMPR7TmUna/m9gd91epPSo/V6/9yBWoIpdJxZeerOfRjBt4hjngHO+o7HEW5Yb0+yStKR2m1On5ZlYyBoYY75vjpOxxFuSm9pp58ws9fUHd6D6KxFo22FgNtHRptHdr65yk49C1HMvYgNQZIYYjW2AqdiQ3CzAaNmS2Glo5YOvXD3tUbexdPNAbqBltfdvC7MxRkVjD5iSFdvmyk0vc01NdSlJ9F6blsaopz0VaeR1dVArWluN/5DP38Qjv1er0myded2UdgwWbqMaZeGNMgTGgQxggkJo1V2NdkoEGHoWzETFZjJaswFLqrztMgDSgQ9hQbu3LBagDSwRcztwCcvINx7eer3gB6uZT9+ST8ks2QOzwYpFZ+Um6QtrGB/DNJFGWdoCY/BU1xOjbVWTg15GEny3ETkitv5V+QZmTk3wWdnOR7xELeN2PF878SOMad2x7wvWy71OmoulDOhfJiKovzqSrMoa4kB135WYwu5GFZcxa3hixsqGo5plqakGU8iHK7IRh6DsPJbySeg0IwMOw175V9WnZyMT8uP4bbQFum/WmoutmqtEtjfS05KfEUpcWhy0vEtjyZ/g0ZmIqGln2KsKXQ2IsK837orNwxtHXHzMELC0cvLB09sLFzxtjE5KbiuNZC3n02OwmNBktrOyyt7XD1ar0MsdTpKCnK53zGCSpzk5HnTmBTmsTQ85swK1gHR5refc+YBVHtOhKbwXcwYOjtmJiqAlY9TU5KCVs+Po6diwVT5g1RCV65pqqKUs4k7KQyNRabgsMMrE/BRzTgA1RiRo6JL4kO9yPcgrH1CsJ9YDCOtg7o6/Z9r0/yAm54eI3QaLB39sDe2QNG/ndhq8aGejLSjlKUehBdziGcSxMIyfwIMj+ibosRySYBlLuMxC7kLgaFjsXQyLhTnovSNU4nFLBtZTI2TmbMeD4UUwu1wIzyX7W1NaQe+oXqk7/gWLAf74bTBAsdWinINBpIostMDPqNxMV/JJ4DAgnsZl26vT7JdwVDI2N8AkfgEziiZVtZ0TkyE3ZQmx6LQ/FhIrI/xyDnMyp+tCDdMpxGn/H0HzkNF8+BeoxcuZSUkoSt2ezfdBoXH2um/jEEMyv1htznSUl+RhI5cT9gnLUL3+oEQkQdDdKAVOMA4jznYuk3loHDxjHQ2o7u/j9aJflOYuvoSujEh2DiQwCUF58n/eBmtKnb8C47gPPxX+H4G2RovDnvPgHHEfcxMHg0QqO6BfShprKeHf9OIfNYEYPCnZnwu8EYGnevFphy62i1Wk4e2k5lwia8CnbiKfNxA3KFK0lOd2M6eBIDR0whyLrnLfvY+5O8EEg9TIeycXBh+NTHYerjSJ2OjJR4zh/ZjFX2dkbkrMQg93PyNzmR5TwOq6H34h8xUXXr3AJSSjKOFvHrV6eorW7gtgd8CRnnqapL9kH1tTWcOvAj1ce+Z2BJLEMoo14acNJsGDn95+IxYhpeA4PwFD37b6OjpYZXAvcABVLKIc3b7IF1gDdNpYYflFKWNj/2MvA4oAX+JKX8udMi70GERnNJ985blBSc5fSeDRilbWHY+U2Y/LKe0l+sSLMbi3nYLAaPuluN2OkC5YXV7F6XRtaJYhw8LJj2p1AcPS31HZZyC9XW1nBy9ybk8Q34VewjmBqqpCkpVpFkD56G/233MdTGXt9hdqqOZpJoYBnw5SXbFgDbpZSLhBALmn9/SQgRCMwGggB3YJsQwk/VlKfpZu59zwHPUVVZRtKeTehObiaoZDsW23+kaLst6U4TsY2Yjf/w8apL5yZdKK3jyE+ZJO3Nw8BAw5j7BxE8zlONoOkjtI2NJO/fQlX8OgaX7mSYqKIMS07YTsAkeDqDR9/DcDMLfYfZZTpaTz5WCOF9xeYZNC0kAvAFsAt4qXn711LKOiBDCJFO09KA+2883I4T3bx4jYWVLWF3/R7u+j01VZXE/7oBkbSRYQXfYfLjN+RtcSbLbQouox9hwJCR+g63R6koquHojhySYvOQOknAGDdGTPXB0u7mxiMr3Z/U6TidGEvxgTUMKPiFYEqbWuy2YzEOfZCAMdOJNO4bfwed0SfgIqXMB5BS5gshLk4T9AAOXLJfbvM25RrMLKwYPvX3MPX3VJQVc2zXWkxSNjHi7GoMN3zJ6U0DKBz4G/zu/H3TsE7lKlJKzqaWcWxHDhnHitAIgV+kKyOmemPtaKbv8JQulpeVRuaOlXhlf8cgmUc/acgJi0iyh9xPUNT9DDe30neIt1xXdvy2drei1Ta1EGIeMA+gX79+XRhSz2Ft68CIe58FnqWk4CxpO77EPn0jkal/o+HUEhIsIiH0IYZEPYBRH2mRtKWypJbUuHOcOnCO0nPVmFoaMXxKf4aM9VQt916utvoCx7d/hcmJtQypTcBdSJKMgzng/xQB4x8mzK5vVxHtjCR/Xgjh1tyKdwMKmrfnApcube8J5LV2AinlCmAFNJU16ISYLj9/Z5/wFrN39mDk7JeBl8lIPsT52FUMOvcjjvueoWTfy6S6TMX59sf6XHdObVUDGUeLOHXwHGdTS0GCu68twyb1wzfcRQ2J7MWaumN2UbI3moDiXxhBNfnCiUP9/0D/CX8gqH+AvkPsNjojyccAc4FFzd+/v2T7V0KIJTTdePUF4jrhen3axVE6jfV1JO7ehO7IGsLOfYPxhq9J3zSQIv/ZBE7+A9a9bITARZUltWQcLeRMYhF5aWVIncTayYyIe3zwH+mqumR6udLCPNJ+/heuZzYwSJdLjTTmhG0UZhG/IzByKm7dbLZpd9DRIZRrabrJ6iiEyAXeoCm5rxdCPA5kAw8ASCmThBDrgWSgEXhGjazpPIbGJoROmA0TZlNWdI4jv6zEOX0dkcnvUJW0mDj7iQROfw5L7xHNd597Jp1WR0FWJdlJxWSdKKYgqxIAO1dzwib1wyfUCef+Voge/ByV65CS0uQdnPlpGcEVsUSIRpINAzkQ+CaBE+cyopc2aDpLr69C+dlfY/Eb4crY2b1/8YeWEQWx/yK4dDvmoo5cU19qQ36Hz7i5GJjZ6DvEdqkoriEnuYSc5BJyT5VSV90IAly8rfEZ6siAUCfsXHvvkDelSV1FIRnbPsP25BpcG3Iol+acdL4b53FPMSDwqmKLfZ6qQtkHCI2GQWFRDAqL4khqJlk7ownM+xb/uNeojnuHdOfJuIx/Chf/Ud2qdV9dUU9eWhl5qaXkpJRSdr4aAEs7EwaEOuEVaI9XgD2mlqpwWK8nJXnHdlD867/wK9lBAA0cxY/9/V/Fb/xviezvqu8IexyV5HupMD9vwvzepLb+Vfbu24bu0ErCz/8Hs6+/J9PYl6qhj+E/YS6Gpre+RVxZUtuU1NPLyEsta0nqhsYa3H1tGTLWA6/B9ti5matumD6i/kIpqVtXYJu8Gs/GbKykGfttp2Ix+gmGhY9hqJq4dsN6fZIXCOhmXVK3kqmxIWOipkDUFM6eO0fKz5/hnfE1QYdepvzQO6R5zKTf5Gdx7tc1oxGklFQU1TS31JsSe0VRLQDGZoa4DbJh8Bg33H1tcepn1WNmoVZXlFKYf5q6smK0GrC2dsKhnx+mZqpMQkfkpcZzbtuHDC7YwhDqSBa+7PB7nSGTHyXKwUHf4fUKvT7JK//l4eqKx9xXaWh8mbjYzXDoU8Jy16D5fDXHLCIxGv0UAaOnITQ3PkJBSknpueqmpN78VVVWB4CppRHuvraEjPPC3dcWB09LND2gMJhWpyXl0E+c3fYjjUdPYJ1dgkPZ5WMIKoAyAWW2hpT5OGIYFsKgCTPxGToWjSpLcRltYwNJO9diePgzAuuO4iCNOGw1AZMxTzFsZBSBPeBvoifp9TdeP//rbnzDnRk7x7/TztmbnM1MI3PrMgbnfYs9FeRqPDgf8FsGT3kSc+vrj1rQ6STFZy+Ql1ZGfnMXTE1l07Jn5jbGePja4u5ri7uvXY/qfmnUNnJsx3rOfbce+7g0bCqb1gMucDSiyscZBvTD3NUDA2trNFJQU1FCw9mzaDOzsU8vxK686U2g0MGQipGDGXDvQwTcNr1PJ/yK4nOkbFlGvzNf4yoLyceR096z8Z3yNC6uagb3zbrWjdfen+T/ZzeDhjtzh0rybaqprubo1mhsjkczWHuKakw46XQ37hP/hJvfsJb9tFodhVmVLX3q+enl1Nc0AmDlYIqHry1uzYndxsmsxyT1i7KTDnDi38uw3pWIQ5mWekPICXbGcuxY/Cc/iNuA4PadJzWeUz+tQ7tzLx6nSjDUQYGTERfuHMGwuX/G3XtIFz+T7iP7xD4Kt3/IkJJfMBENHDceSt2wxxl65xyMVHntTqOSvEry7SKlJPnwr1Ts/piw8u0YIIjjLkoc70fb6Mm5jAoa65tatXau5k0JfVBTUreyN9Vz9DemtvYCB9Yvo3bDd/RPLUcnIDvADoupUxj2wFNY2Tpf/yRtKC3MJX79R+i2bMfrdCU6ATmDHbCdPYvwmU/2ynUEtA11JG9fjcmRz/CrT6ZampBoPwWn8c/iGxyh7/B6JZXkVZJvl/raRvJPl5OXVkZWUgEluVVIqQF02BiexcTFkMDxo/EJ9sDcumcnp4KMZOI/XYTtL/HYVuoosTWgfHIEQ3//Fzy6qKWdlXyAE19+iN2OROwqdBTbGlB+10hG/GEBzh6+XXLNW6mi6BypP/4Tn4yvcaCUHOFG5oCHCLrrKewdb+7NUmlbnx0n391LDetbY4OWc2cqOHuqlNyUEs5nViJ1Eo1G4NTfitCJ3jh5m1OY9RP2SSsIaDxJ5S9mHDt2L/2n/hm3/j3rzVNKSdr2TWR99hHuR/PoJyEnyAHDhx4mcsYfMDDs2rH4/QMj6b8okvr6Gg5+8xHV6zYyYO0+8tdPJy7CC5/fP03Q2Hu7NIaukJGSQMEv/2Bo0RbCRQMJxuGcCX+CYePvx0stgKNXvb4lv/KF3Qwc5swdD/WsZNRVdDpJYVYluadKyE0pJf90OdoGHUIjcO5vhWeAHR5+drgOsMHI5OpRNqcO76Dy1w8JrdiFQHLM6nbMxs4nYMTEbjXB6kraujqOrPk/atasx+lsFRXmkDc+iKGPv8CAwfot7JaZuJukz/6OW+wpzOohz8sco4d+w6iH/oyRSfetxSN1Oo7u2Yzct4xhtQepk0Yk2k/G4c7nGRQ04vonUDpVn+2uUUm+afGMrBPF5Jws4WxqWcuNUnt3CzwD7PAMsMfd1xYTs/a3uPJzTpOx5Z8E5W/EhirSDQdRMfQJQiY/iqFx9+mbry7M59DH72D6wy6sK7XkORtSdd94bv/9K9jZuOg7vMtUlJxj/6pFmG7agXNRA6XWBpRPH83Iea9i79x9SnDX1taQ+J+VOB7/lEG6DEqwJr3/HHzvfg47tc6B3vTpJD9gmDNRfSjJaxt0zX3qxWQnFVN6rmlGqZWDKV4BdngE2OHpb98pfepVleUc2/Iv3FO+oL/MpRA7MnzmEHDPc1g76G8K+rkTh0ha/j6OsUkYN0KavyXmj8xm7Mz5GBt273sJjY0NxG36hMov19AvrZxaI8gZ60vQH19k4JDb9BZXceE5Tm5eil/WWpwpIVvjRVHwHxhy1zyMTc31FpfSRCX5Xp7kK0tqyTpeRFZSCbkpJTTW69AYCjz87Ogf5EC/IHtsXbpunLpOqyVx17do4j4mtC6eWmnECccpuE76C57+YV1yzStJKUn9+RvOfrYctxPnqTeEtJEeeD/xLMNHzuhxwzkBUg/+TPq//onXgUw0OsgMccT18XkMm/jwLRtzn5F6jPyf/0Fo0Y+YizqSTMMQo59h8G0zb2rinNK5VJLvZUleSklJfhUZiU211Quzm0rwWjua0i/Igf5BDnj427Xar97V0pMOUbBtKcNKfsJM1HPcbASaMfMJHD2tSxYl19bUkPDvf1KzdgOO+dWUWQrOTgpm+JOv4NN/aKdfTx8Ksk8R//FCHH+Kx7JGctbTDKM5Mxn1yF8xNun8VrTU6Ti2/2ca9nxIWPU+GtFwwn4SjhP/Qr9ANQSyO+ryJC+EyAQqAS3QKKUMF0LYA+sAbyATeFBKWdrWeTo9yb+4hwFDHYl6uOevFCN1kvOZFZxJKOTM0ULKC2oAcB1gjc9QJ3yGOnZpa72jCs+f5dTm/yMgZx2OlHHawIeS4HkMvesxjE1uvt++6nwe8R8vxOyHWCyrtOS4GlJz/52M/d3L2Fn3zuF61ZWl7Fv1Hkbf/IRzYQNl1gZUTr+dkU+9io3jzfeH19XXceSnL7A/ugJ/bRplWJLm9SCD7vkzdi7d576AcrVbleTDpZRFl2z7ACiRUi4SQiwA7KSUL7V1HpXkLydlU2JPO3Se0/EFVJXXo9EIPALsGBDalNgtbLr3Gqa1NVUc3fIZLkmf4q3LoQB70n0eJvCe57B1cOrw+c6dOETSR4tw2p2MQSOkBlph9duHGDvt6W7f395ZtNpG4jZ9QvmXq+mf2tRvnx8VSMgzL+MZ0PFa6yXFhSRtXsagjNW4UUSuxp2CoMcJumseJubWXfAMlM6mryR/Coi6ZP3XXVLKNvtNVJJvUpx3gbRD50k7dJ6Kolo0hoL+QQ4MDHPGO9gBE/OeV1td6nSciP0Wzf5lBNUlUCVNSHSajuddf6H/wMC2j5WS1F82kvvpR7gfP0edIaSN8sT7ifmEj5jWbT696EPSgR85869/4n0wF40OckJd8Zr3LIHj7rvu65J5+iS5//kHoYUxWIoaUkxCkJHPEHDHA6q/vYe5FUk+AyilaerRv6SUK4QQZVJK20v2KZVS2rV1ns5O8qte3IP3UEfG9YAkX1Vex6kD50iNO0fx2SqEAM8AO3xHuDAg1KlHJvZryTxxgOJtSwgp3YYGHUcsbsdk7HMEj5xwWWLS1dcT/9VSqld/jXNuFeUWgtxJwYQ99TID+ofq8Rl0P2czjnPk43dw3XoMy1pJXj8LzH83h4hZ8zG4pHSClJJjB7dTF7uU4VWxSAQn7CZgP+F5+gXrb/SOcnNuRZJ3l1LmCSGcgV+A+UBMe5K8EGIeMA+gX79+w7OysjolJuj+SV6r1ZF1vJiTe/PISipB6iQuPtb4RbgwaLhLjy8dcD0l+Rmc2bwE/7MbsKKaJINAKoY9SWDEZBI+fx/TTduxqWgk38mQC7/pnuPbu5uKskL2fr4Qi2934FTcSKmtIbX3TSBs7kskH/oJm4R/EaQ9SSXmnPL4DQPv/gt27gP0HbZyk27p6BohxJvABeAJ9Nxd012TfEl+FSf35XPqQD41lQ2Y2xgTEOnG4NFu2Lr0vTHHtRfKSN6yHNP4VeSekTilGGHaAOmDzDH/7Wxuv+9PGBt173sP3U19Qy371i+levU6fDKqqTGGwoAGvHw01A3/PYOn/hFTS9vrn0jpEbq0do0QwgLQSCkrm3+eBLwNxABzgUXN37/vjOt1LDi6Te0anVZHxtEiju/K5WxqGRqNwDvEkcGj3egXZI+mh6yK1BWyE/dS8MNOvA4b4ibglJ8Gb98SbjOvILE4m7yzZ/H2Vq3NjijMy8WwIpMxEXmkB0Fmmj3+x43QHYeCvJ1YuA/F77ap+g5T6WKdVTnIBdjU3JdqCHwlpfxJCHEIWC+EeBzIBh7opOv1KDWV9STtySMp9iwXSuuwsjdl1MyBBIxy6/XdMW3RNTZyaOPHVH6xGo8zFTiaQNqUwYT88WXu9w0nO3EHZ3f8gzvO/5uGVV+x22oCVlHPMXT46D59o/V6kuN3U7nzn4RV7sQVHcdt7sD+3ucZPmwcZ1LjOLZ8EV47T6L9w1/ZPuAtrH/3CMPv/yMaVUisV+r1k6FWvbQH72BHxj1y67trinIrSdyWQ9rh8+gaJZ4BdgRHeeId4tgjlr3rKtWVpRz4/D2MNvyEY1FDU7nde29n1BOvYu9w9VjvkpyT5GxZjH9+DKbUc9gwjJrwpxk54T6MjdQIEACtVkvitq8xPvwJwQ3HqMKUZNd78b77Lzh5Xd1DWlycy/5PF2Lz3R4cy7QU2xvReP9kRv7hFcys2xwboXRTfXbGa/RLe+g/xIFxvx3caedsi5SSs6dKSdiaTXZyCUYmBgREujIkyhN7N4tbEkN3VZCbSvxHC3H8+TCW1ZJcLzMMH76PMXP+inE7qi3WlheS9uM/8Uhbjb0sI5X+ZPn/nvC7n8DOum8uoF1dVcHRzZ/gkbKKfjKP88KB7EG/I/CeP2Fh8//bO+/wqsr7gX++udmbDJJASEhCQAgjDNk42CKgddSNE/do1bba/nxq/alPrbZqW6vitloQ/SEqiDJUZMoSMKwQRiAhCUlIQva49/39cU5qjAFyQ+69J8n7eZ7z5J5x7/nkjO95zzvPPHxjbW0Va+f/jYb/LCLhSDWV/kLxtGEMu/t/iE60VjmW5vToIO/iIO9wKA5+X8j3y7M5nl1OQKgvQybGkzahJ/5BnafqY1vYu2UF+1/5G4kbDmOzw+HB0cTeOrfN/a+o+moyV75F0LZXiK/P5rjqxo64K+gz436SErpGq8yiY9lkLn2e/rkf0Y1ysrz7UDHsTgZNufEn1SVbi1KK0J9yCQAAG29JREFULaveJ++NeaRsL8QhkDMqkdQ7H6LPqCku+A807Y0O8i4K8na7g30b89n2RTZlhdWERQeQPiWBc8bE4t2FsxLsDjvfffY6JW+/Q/KeEmq9Iee8VPrf9RtSBk1on50oxdEtS6he/SJ9KzZTo3zYGDKVsAvvJ33YqE6Zb38oYwNFK55nSOlKvHGwM3gsAec9QL+RU9utP/8Du9eT8fIz9FqdSUAd5KSG0+3mGxl2yVy8bF33mrY6XTfIP7KOhLQIJrZzkHfYHWRuKmDz0kOcLKohOiGEYdMSSR4a3aXz2ysrSlj39p/x/ugL4vLrOBnkRenMMZx712NExCa6bL8lh3eQs+yv9C34HD/q2ew9nJoRdzJq0uUdPt9eOezsWv0RXhtfYkDtDiqVHxnRM4m/6EF6prhuQPCi49lsnPck3T5dT8RJB4XRvjh+eTGjb/4d/sFhLtuvpm3oIN9OQd7hUOzfbAT3suPVRPUKZuSsZHoPiuyUKcfWcuzwLra++jTdv/ye0CpFfpw/tqsuYdSND+MX4L788prSfPZ//iLx+9+nmyojiwSyU+cwbOYddAvrWH2wVFecJGPZq8TueYtejlwKiORA0nWkzbqPsAj3dcBWU1PB2veeQ83/hPjcGioChRPTRzLi7seIjE9xm4fm9Oggf5ZBXjkUWVuPs2nJIUoLqoiMD2bkzCSShkR16eCese5TDr3+jx/7TRncnbhb5jJo6rVu6++8JVR9DZmr3iZg66sk1B+kWIWyPfZyki96gKTeSR7zag35OQc59PkLDDj2EWFUkmlL5cSQ2xk6bQ5+7dB7Z1txOBxs+fIdjr/1Bkk7i7HbIGdMMv3ueoiU4RM95qUx6LJB/p1H19GrfwQT57Q9yOfsPcH6RQcoPFJORI8gRs5MIjk9Gumi2TIN9XVs+PDvVL2/kIQD5dT4wLEL+jPozt+SkDba03o/RSmObvuCim9epH/5BmqVN5tCJhNy/v0MGWGd+vbK4WDvllVUrHmZ9JPfYMPB9uDx+I2/lwGjprqkH/6zYf+O1ex55VkS1hzArwGO9I8g4qabGD7rVo8+3LsyOsi3IcgX51awftEBjuwqJribH6MvSabvyNguG9zLTuSx8bWnCFq8msiSBk6E26i89AJG3/YHQqPiPK13RkqO7OLI53+lX/5n+FPHTttAqobcxNBp1+PnoQGza6rK2fXlm4RnvE2K/SDlBLA3ZhY9p/2aHsmn75nTChzPO8DmV54kcukmwioc5Mf4Un/ldMbd9AhBwbq+vTvRQd6JIF9RUsN3nx1i74Y8/AK8GT69N4Mu7Nlla8tkZaxl96vPEb96HwF1cDQpmMDrrmL0L+/D27fj9SdTU1ZI5rKXiN73H+JUAcWEcSD+FyROuZuYRPeMIHZ4fwZ5K//JgIJPCaOSg16JHD9nDoNn3EZgB+xPpqa6nA3/fg7mf0psXg0nA4X8yYNIv/0REvsM9bRel6BLB/n4/hFMakWQr6tpYNsX2WxfdRSlFIMv7MXw6Yldsp57g72BjZ++Qsn780nOOIHDC46e24vEufeQNv4ST+u1C8phJ+Pbj6n/7nWGVG1EgJ0B59KQfiODLry83VP3NZUn2b3qPfwy5pNWt5N6ZWNnyAR8Rt/BwDHTO0XfRQ6Hgx++eJ/8t98gfmcBSuBQendi5tzMiGlzdFaOC9FB/jRBXilF5qYCNizKorKsjr4jYxg1O5nQKM+8wnuSE0W5bHjzaYI+/ZaYogbKg7wonXouQ+96lOiEzjVOblOOZe/nyIqXSclZRDQllBJMVtRkAkdcQ79zp2BrY/3whrpa9mxcRtW2haSVfEWwVJMjseQm/oKUaXcSFde7ff8RC5GXuZ0drz5D1KodBNUojsX50nDZVMbe+AghoZGe1ut06CB/iiB/PPskaz7YT/7BMronhjDhqr7EJne9OsB7tq4g8/UX6bXuAAF1cCwxGP+rLuPcax/A17/rdH1sr69j15qPqd22gLTydQRKLQVEcDRiLF6pk0kYPp2o7qcvfyg6dpDD21bi2L+KfmVrCaOCKuXHrm4T8Rsxp9Ok2ltLTUUZ373zLOqjpcTk1VDhL+RPTGPg3IdJ6j/K03qdhq4b5H+/jvh+3Zh0408LsapO1vHdJwfYvT6PgGAfo1fI0XFdqlC1traKDQv/TvUHi+idVU6dDY6NTiLllnvpO053QVtRXkrm6oWw9zP6lG8hVKoAKJAo8vySqQ+Ipt43DBsOVG0FwdW5RNUeJRZjBMyTBLEndDw+abMYMOFS/ANDPPnveByHw0HGqoXkvjWPXt/nIQoODYok7NprGD37dnxsXS9btD3RQd4M8na7gx++zmHz0sM01NoZPDGeERcn4RfQdbpZPXx4B9vfeI7oL7cRcdJh1JKZOYERtz3i0lapHZm6ujoO7VhNyZ5vsRXuIrLqIMH2MkJVOXZs1Ig/Rd7dqQhMoD52KN36n0+fQaOxeevA1RL5BzP4/rVniFy+jZBKB/lRNkpnjmXkTQ/TM7avp/U6JB4L8iIyHXgRsAGvK6X+fLrtXRnkj+wuZu3C/ZTkV5GQFsH4K1PpFts1eoasa6hj3bLXKFmwgD7bi/CxQ845kYRfew3DLrtdByONR6irrmTrf/5O7QcfE3PEaHNxYEQcsdfcwKhJ1+vUvRN4JMiLiA3IBKYAOcBm4Bql1O5Tfae9g/y7v19PWPcAvH1tHN5ZRGh0ABOuTCWxC3RDoJRi96Hv2P3ev4hcvo24IjvV/sKJC4aQNvch4tJ+dj1oNB7jyKav2PfWP4leuxe/esXRWG8qZ4xjyPX3kdwjzdN6lselw/+dhpFAllLqoCmxALgEOGWQdwU5e0vw9rMx5hcpDJnYC5tP5y70arA3sPj/nqZ20RIG/lDOQDvkJ4dTfvOlDL36HnyCumbf6xprkzByIgkjJ1J7soRt776I/8dL6PXmasr+vZr3hkUQdtUvmXnR/Z0+cdbeuDolfwUwXSl1mzl/AzBKKXVvs+1uB24HSEhIGJ6dne30vg7v2EbBwayfL/+hGHuDg179I/D1b1YNzsmLxR0Xl9P7aLJ9dXkJOVtWY9t9kODKBhq8vahLjSd6xDiC4noZm7dNysnNnd2Ls7/v5M+36UsudnL2mLblzDl93pzegZO/79z2SinKDmdSuH41AVm52OyK0m6+kN6f3qOnEBDYUodzLr5W23IanPhS4uChhERGOb8TPJeSb+m/+9lTRSk1D5gHRnZNW3Z0YOt3bP9y6SnX5+1ry692UMLCobEWaFkVR1at8KiORnPW9Ij+8XNWHjlZ73rOxYVc9sjjbQ7yp8LVQT4H6NVkPh445oodXTBnLuffcFvrv+DkG4z6+bOp3ffhzC4OZqwjc8HrhK/5gfByB1V+QsGYFJKuvY1+I04xkk8b3tqc/4qTx9WFx+jHrzi7Dxf/D87SpvPmWifnz1sb/ocWzptyKDI3LCN38QfEfJ+NXx2c6OZN5YR0Blx+CzG9W9/fj/PHqE0Xn1MEhLV/Gx1XZ9d4YxS8TgJyMQper1VK7TrVd9q74LUzUVacx+b3X0B9vor4w5U4BI6kRRJ8yWxGXnG3W/tt12g8TeXJE2xe8HdqP1tGwv6TOICcft0IuPRiRlx5T4fsA+hs8GQVyhnACxhVKN9USj11uu11kP8pddVVbFv6FkWffkyvbbn4NkBBd19qpo1h6PUPEJPongHKNRork71nEz+890/CV20jstROlR/kjuxN7KVXMmzq9Xi3YdzbjkaXbQzVEWmor2PHivnkLV5IzOZDBFcrKv2F/LF96H3NLQwYN1t39KTRtEBDQz3bV7zP8Q8XELc5G/96KA3xonjsOSRdfgP9x3fee0cHeYvjcDjYteYTshe9R+T6vYSXO6jxgdyhPYmcdQnDZ96Kb0DX6UNGozlbqitK2bLoVcqXfk58xnF87FAU4UPF+UNIvuJG+g6b1KmqY+ogb0Ea6uv44esPyf38Y8K/20dkSQP1NsgZ2J3gi6Yz7Be3Exyme+vTaM6WksIctn74L+q+/IrEzDK8FOR396FyzEB6z7qKAWNndfgUvg7yFqG66iTbl75D8fJlRG87TGilot4Guf0j8Zt0AcOuvIvwqJ6e1tRoOi3Hc/ez84NXsH+9jp5ZZdgUlITaKBmZStyMSxk0+Sp8fD03lm5b0UHegxzN3Mq+ZQuoW7+JuN3H8a+Hal/IG9KDkCmTGTr7VkLCu3taU6PpcpQUHGH74tepWvUN8bsL8W2AigAhf1AcQePHMXDG9XSP7xgdpukg70Yqy0+QsXIhhV8vJ+T7LLoX1gNwIsxGSXpvoqfPJH369brKo0ZjIcrLiti+5C3KV64iaudRQiodAOT18Kd2RH9iLpzGgAsux9+i960O8i6ktCiXvd8spnjjt/hmHCD2SCXeDqjzhmOp3fAaPZyUqZeTPOS8Dp/vp9F0BRx2O/s2fcHh5YuRTTvoeagcbwfU+EBB71Acg/sRM+ZC+p93KYGh1hiwXAf5dqK2uoKsbV+Rt3UdNbt2EZiZS0xeDV5AgxfkJQRRPyiV6HEXMHDyVV2uQYZG0xkpPXGMPSs+pGTtt/jvPkjMsRq8FNgF8uMDqTqnF8EDB9Nr+PkkDR7vkQHudZB3krraKo7u20r+7i2U799Dw6EjBGQX0D2vBm/jLY5Kf6EoMQzH4H7Ejp3IORNm66Cu0XQBSouPsWf1Yoo2rsFv1wG6HynHz8iVpdYbCnsGUZMUi29SEuF90+g5YCQ9Ugbj5e26nmQ6fZA/dvAHKkuO4xsYgn9gCH5BofgFhuAlXtTX1WJvqMNhr6e+tpqK0kIqiwuoLi2ipvQENcXHqc/PQx0vxqeojMCSasJLG/4bzAHKgr0o7RGCvW8iIYPTSTx3Er36jdDZLxqNBntDPYd+WMfRLd9QkbET36wcInMrCKr5Mb7W2aAk0ofqbkE4osKxdY/GJyYG35hYAiNjCImIpUfq0DZXm+70QX7J3bNJ+Wp/m/fb4AVlYTaqwgOojwqF2O4EpvYlql86CQNHExZ5+sGbNRqNpilKKYrzDnIkYyPF+3ZSd+gQXseO41NcTlBpDaEVDryahd+yJ+9l9BX3tGl/nupq2G2kzLmLwtHbsVdXG1NtDY7aGgDEZkNs3oi3N2Lzxic0DP/wCALCowmKiCY0uieRccnYbJ3mcGg0Gg8jIkT1SCGqRwpMve5n6+tqqyk+doCyvGwqThRQdaKAc0ZNbn+PzpKS12g0mq7MqVLyOkNZo9FoOjFnHeRF5HERyRWR7eY0o8m6R0UkS0T2ici0s92XRqPRaJyjvTKhn1dKPdd0gYgMAK4G0oAewEoR6auUsrfTPjUajUZzBlyZXXMJsEApVauUOgRkASNduD+NRqPRNKO9gvy9IrJTRN4UkcY2vj2Bo022yTGXaTQajcZNtCrIi8hKEcloYboEeBlIAdKBPOCvjV9r4adarMojIreLyBYR2VJYWNiGf0Oj0Wg0LdGqPHmlVKsqb4rIa8ASczYH6NVkdTxw7BS/Pw+YB0YVytbsS6PRaDRn5qzryYtInFIqz/z8a2CUUupqEUkD/oORD98DWAWknqngVUQKgew26kQBRW38rivRXs6hvZxDezmHVb3g7NwSlVLRzRe2R+2av4hIOkZWzGHgDgCl1C4RWQjsBhqAe1pTs6YlydYiIltaagzgabSXc2gv59BezmFVL3CN21kHeaXUDadZ9xTw1NnuQ6PRaDRtQ7d41Wg0mk5MZwvy8zwtcAq0l3NoL+fQXs5hVS9wgZvlOijTaDQaTfvR2VLymlMgIi21W9B0MPR51DhLhwryIhJnxYtcRHqIiPsHdTwDIjJIRH4HoCz0yiYisZ52aAkRifG0Q0uISD8RuQgsdx4TRSTB0x4tISL+nnZoCU/EsA4R5EXET0ReBlYD80TkMk87AYhIsIj8DVgGvC4i15rLPXpcxeA5jHYK3iLi40mfRkQkQEReAL4QkefNFtMexzyPzwPLRORVi11ffwXmA76e9mnEPI/PY1z374jIXeZyj8cTEQkSkXnAH0Uk0lzm8YShJ2OYx09KK5kNxCml+mK0qH1CRPp6UkhEegBvY9x844BPgMZUs+PU33QL0UAcMFwp9ZRSqt7DPo3cA0QrpdKBxcDTItLHk0Ii0hP4N8a9MAPjJvyLJ50ARCQUWASMV0oNU0p94mmnJtwP9FBKDQAeB34Fnr/uzdT7E8B4IAS40PSywtuPx2KYZYO8iAQ3mVVAIYB5sX8B3CEi4R7wCjE/lgEPKaXuVUpVADHAYhGJNrdz67Ft4gUQhtG6uE5EponIw57qz7/xPIqIDeiG2e2FUmo1UImR4grzhJtJDfC6UuoBpVQ+sBDYLiKDPejU6PVvYBeAiIwTkakikmrOu/3eFRGbuV8BdpqLewBLReQcd/s08Qo0P9Zi9KV1HrAfGC4iKeY2bk/NWyWGWS7Ii0gfs6Xs2yJysYgEAdXASTP1DPAsMAyjr3q3nMDmXoCPUipbRAJF5AHgESAI44IfoJRyuNnrLfN4RQAVwDoReQL4LUbAeEFEbmx24bnD6x0RmWkuLgdGicgQ82G4F+gLJJvfccfx6icir4hIAIBSqhj4pskmvUyffa52OYNXHfAVoEQkH3gamAKsFpE0N15f//VSStnN1PoxIEFE1gDPYJzXlSIyxZ3BVERSReRdjOyP2UCIUipLKVUEfA3444HUvOVimFLKMhPGQ2cJ8BhGf/QvA38G/IClwEWAr7nt48BHHvJ6CfiHuU6Avk22fQJY4SGvfwHPmev+gRG8hpjzVwAfYdwI7vZ6BaPls4/590Ngu3k+/xeY56bjNR7YBDiAPzSev2bb9AMWucPndF5NjuMk4OFm19cXFvAKw3jriTWX3QN87sZjdgNGlyl3AbcArwFzmm0zF3geI9vSXV6Wi2Fuu5BbeYB6Au8Btibzm4BRwJXAW8BIc9055on18ZDXBmC2OS+NwQIjFbgYCPCQ13cYr6tDgBXALU22/xojL9VTXlPN+SQgwvx8OfDrxuPoYq/+wECgD8YgNoktbHM18Kz5eS4w2A3Hq7lX7ybr/Jttm4qRV+/vKS/zeu9pBtBkc5kfRqIi0tVe5v6mArOazD8D3Gl+9jb/JgD/A9yN8aZ9nhu8LBfDLJVdo5TKBUZgvJY2zv8L+JNS6kMgE3hURB4CFgAHlRsKFU/h9TLwa3NeKaWUiIwB3gTWK6WqPej1mFJqB0bruVlijLW7BsgATnjI6yXg9+b8IaXUCRE5D3gQc3AZZV75LvTaA2QppbIwHoBPwM/ytycBkSLyf8C1GFldLqUFrz+ZXqKU+u/+RWQs8Aawselyd3uZ5ykf44EzV0RuAr4ENmOUVbkcpdRyYLmINPa/VYNRPoBSqsH8ewQIBp7EeHi75NpvmtViyRjmjqduC0+7OGBss2WNT76bgLVNlodjvBaei5GCmAC8CFxvAa/5pk8Qxg3wPfBLC3h90Lg9Rp7fQ8DVFvCaj5mawkjB7weudYdXk3WNb1whGKnTSc3WL8Mo7LzCKl4Ygep3GFlcV1nIazBGKnmpK66vM7k12+594LJmy87FGMjoOhd53dpsmcdjWIuu7thJCwdol3kzDWt6ITUeKIwCp181WfYOMNDKXkC6Fb0sfLzCPeHV6Gb+/RWwxPx8jXkDXmBBL2+alPtYyMsdWZJncvMCAoGPMWq4CTAN8HOh02PAVuDBU6z32D3Z0tQe/cm3GvO1OBDj9XwncL6I7FNKVYqIl1LKoZSyi8hvgI9E5DBGKmYgRuGPFb0AUEptt5iXVY+XAlBKlbrZS5R5t2EeG6XUCyLyRxEpwyiU/lgp9Y0FvRYrpTIt5rUQ+KTZdm53U0YtozCMguCLgTsxspZWu8JNjEaGd2CUKdU3WyfKwO335OlweZ5803qg5s1fgfHaVwj0BoY3rjO3tymltmKkGoYCtwMPKKV2W9Qrw6JeVj1euzzkpZpsp0QkTESexQgiM5RSt6p2zOduZ692K99pR6+5Sqmq9gyibXEzGQFcAEzHaLvyB6VUTXu5yU/rsr+B8XbhLyJ9ReQWERnd1MtM6Lj8nmw1rnpFACKBfwLLgfuA/ubyPsDb5ucHMUrFfwXENx4jVzlpL+3V5PteuKDmjPbyiFsQcIeLve5v4vU4UA+sB/6IMVzpHH7M1nLpte/s5JKUvPlkW4YxVuGTGIN43w2gjJL6k2K0gOwHPIAxLmyOub7dX/20l/Zq9GpEGSnFnbQj2sv9bmaquVIp9aqLvXpi1MkHo9D0CWCyUupPwL2Nzqa3y679tuCqPPlS4G9KqQUAItINmCRGR1kDgEEY+WwHgOcwWvWlKqX2u8hHe2kv7WU9r7N2U67rL6clr8li9I9TppT638YNlVKficiDGA+obBf5tJl2C/LNCkP2isjRJsvqgT7KKKjYISJrgTVKqeUikoRRh9Ul9X61l/bSXtbxsrJbK7xSVLNyGxEZgNGK+zBGdw/Ww5m8neYTRun31ZyhKhXGq8yfz2Zf2kt7aa+O6WVlt7Z6YSSQL8V4y7jHncfS2anNefJitHLLANIxCj4al3u18DkJo2kvInKZuHDQCO2lvbSXdbys7HY2XhhVNr/F6KLgJVc5tgdtyq4RozP+GcBM1axakPqxCt1g9WNBTRpQJCLLMJoWr227svbSXtqrI3hZ2e0svUowspAKXeHW7jjxWuPd5HMs8Kn5ORm4lR9bpMUA72KcnB4YnQSdxOi86FIXvG5pL+2lvSziZWU3q3q5emrVgcEo1X4Ro8oQGH2AzwdmARsx6q9mYFT6HwXc1+w3bnbFCdNe2kt7WcPLym5W9XLXdKaDIxg9qL0HXAesBO4y1y0EVgFTzPnxwCGadJuJ2W+yC06a9tJe2ssiXlZ2s6qXO6czHaBQjFZdIeb8NIwuY6djNE7YD0xvsv0SYJDLpbWX9tJelvGysptVvdw5nbZ2jVLqJEb9z5vMReuALRgjnhzA6Dt9iojMFZE3MJo+Hzjdb7YH2kt7aS/reFnZzape7qQ1VSg/BtJFJE4ZHQbtwGwYgDEyzHsY/UofUErNUEpVucxWe2kv7WVVLyu7WdXLPbTidScO+AvwaJNla4FxTeZt7n4F0V7aS3tZx8vKblb1ctd0xpS8UioPY8zSi0TkShHpjdGsuL7JNvYz/U57o720l/ayjpeV3azq5TaceBpehJF/tRe419NPJ+2lvbSX9bys7GZVL1dPjeM3tgoxeoZTyhwo1ypoL+fQXs6hvZzHqm5W9XIlTgV5jUaj0XQsXD78n0aj0Wg8hw7yGo1G04nRQV6j0Wg6MTrIazQaTSdGB3mNRqPpxOggr9FoNJ0YHeQ1Go2mE6ODvEaj0XRi/h9tIZ8l4XamlwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "times = pd.date_range(start=datetime.datetime(2014,6,24), end=datetime.datetime(2014,6,25), freq='1Min', tz=tus.tz)\n",
    "\n",
    "solpos = pvlib.solarposition.get_solarposition(times, tus.latitude, tus.longitude)\n",
    "print(solpos.head())\n",
    "solpos.plot();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAD8CAYAAACYebj1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3deXxU9fX4/9e5M1lICFvYAiiLshMSIKLiBiJiXVBs9YNLBdGq2E/BfqtV60LdWq38tPrR6ketez8KiIpbLYpbqbYIggoBZBEE2cKWhIQks5zfH3cyBEhYc2dJzvPxGCZzZ7lnhpk5897OFVXFGGOMAXDiHYAxxpjEYUnBGGNMlCUFY4wxUZYUjDHGRFlSMMYYE2VJwRhjTJRnSUFEnhWRzSKyqMa2ViLygYgsj5y3rHHdrSKyQkSWichIr+IyxhhTNy9bCs8DZ+217RZgtqp2B2ZHLiMifYAxQN/Iff4iIj4PYzPGGFMLz5KCqn4GbNtr8/nAC5G/XwAuqLH9VVWtVNXvgRXAYK9iM8YYU7tYjym0U9UNAJHztpHtHYG1NW63LrLNGGNMDPnjHUCE1LKt1vobInINcA1AZmbmoF69eh32TsMlJVT9sJa0Y49B0tMP+3GMMSYWqlatAschtUuXI3qc+fPnb1HVNrVdF+uksElEclR1g4jkAJsj29cBR9W4XSdgfW0PoKpPAU8BFBQU6Lx58w47mNKPPmLd9b+ky9SpNOnb97AfxxhjYuH7//ovfFnNOPqZp4/ocURkTV3Xxbr76C1gbOTvscDMGtvHiEiaiHQFugNzvQ8n0kAJW1FAY0wSCCtIbR0r9cezloKIvAIMBVqLyDpgMnA/ME1ErgJ+AC4CUNXFIjINKASCwC9VNeRVbFFO9YtrScEYkwRUa3xvecOzpKCql9Rx1fA6bn8fcJ9X8dRGnEhDKRyO5W6NMebwhMNIrUOw9adxr2iONMPUkoIxJgkoCo63X9uNPClEnr71HhljkkEMxhQadVKQ6r45tZaCMSYJxGBMoVEnhWjGte4jY0wyCIcRse4j70ReXDtOtTEmGaiGrfvIU2LrFIwxSUSxgWYvia1TMMYkk3C49qJA9ahRJwVsnYIxJpmo2piCp6LrFKylYIxJfDam4LXqF9cGmo0xycDGFLwVLXNh6xSMMcnAxhQ8ZmUujDHJxMYUPBYtc2HdR8aYxKcatu4jTzk2pmCMSSJW+8hbYmUujDHJxGofecyxMhfGmCQSDu/+MeuRxp0UrMyFMSaJKLp7LNQjlhTAxhSMMcnBxhS8ZesUjDFJxcYUPGbrFIwxycTGFDxmh+M0xiQTtTEFT9nhOI0xyURVbfGap2ydgjEmmVjtI4/ZOgVjTDKx2kces3UKxpgkompTUj0VHcW3MQVjTDKwMQWPOVYl1RiTRGxMwWO2TsEYk0xUdy+69UgjTwq2TsEYkzzcSTE2puAZW6dgjEkqDXVMQUR+LSKLRWSRiLwiIuki0kpEPhCR5ZHzljEIBLDuI2NMkmiIYwoi0hGYCBSoaj/AB4wBbgFmq2p3YHbksrdsoNkYk0zCYcTxebqLeHUf+YEmIuIHMoD1wPnAC5HrXwAu8DwKW6dgjEkS0R6NhtZ9pKo/AlOAH4ANQLGqzgLaqeqGyG02AG29jsXWKRhjkkYoBID4GlhSiIwVnA90BToAmSJy+SHc/xoRmSci84qKio4smOoyFzamYIxJcLtbCg2v++gM4HtVLVLVAPA6MATYJCI5AJHzzbXdWVWfUtUCVS1o06bNkUVS/eKGLCkYYxJcJCk0uJYCbrfRCSKSIW7/zXBgCfAWMDZym7HATK8DqX5xNRzyelfGGHNEtPrHq8cF8fyePnotVPU/IvIa8BUQBBYATwFNgWkichVu4rjI82CqB2xsoNkYk+iqf7x63FKIeVIAUNXJwOS9NlfithpiJrpc3FoKxpgEp9UDzQ1wTCFx+NwXV21MwRiT6KrXUzXAMYWEYS0FY0zSiLYULCl4y+ezloIxJuE15CmpCUUcx1oKxpjE14CnpCYWx7HFa8aYhBerKamWFHw+W7xmjEl8MZqS2uiTgjiOLV4zxiS86JRUn40peMtaCsaYZFA9JdVmH3nLWgrGmKRgU1JjxFoKxpgkYFNSY0RE7HgKxpjE11CPp5BwbPGaMSYJaHXhTmspeEscJ5qBjTEmYVWPfTri6W4afVLA57PFa8aYhGdTUmPEWgrGmKQQHWi2MQVvWUvBGJMMqmsfWVLwmCO7M7AxxiSo6IQYG2j2ljg+W7xmjEl8YZuSGhu2eM0YkwRsSmqMWJkLY0xSqG4p2JRUj1lLwRiTBKqnpGJTUr1lR14zxiQFm5IaI46zu6/OGGMSlEYPx2ktBU/Z4jVjTFKww3HGiC1eM8YkA5uSGhvis5aCMSbx2eK1WHGspWCMSQLVLQW/JQVvWUvBGJMENGhVUmNCfP7d83+NMSZBaTDo/uH3e7ofSwp+PxoMxDsMY4zZLw25ScFaCh4Tvx8CwXiHYYwx+1fdfdQQWwoi0kJEXhORpSKyREROFJFWIvKBiCyPnLeMSSwp/t3NMmOMSVDRbu6GmBSAR4D3VbUXkAcsAW4BZqtqd2B25LL3/JYUjDFJoKF2H4lIM+BU4K8AqlqlqjuA84EXIjd7AbggJvH4U9CAjSkYYxJbQ5591A0oAp4TkQUi8oyIZALtVHUDQOS8bW13FpFrRGSeiMwrKio64mDEWgrGmCRQPdDcELuP/MBA4AlVHQCUcQhdRar6lKoWqGpBmzZtjjgYSwrGmGSgwSCIxP8YzSJyjIikRf4eKiITRaTFEexzHbBOVf8TufwabpLYJCI5kf3kAJuPYB8HTVL8YN1HxphEFwx53kqAg2spzABCInIs7jhAV+D/DneHqroRWCsiPSObhgOFwFvA2Mi2scDMw93HIUlJsZaCMSbhaSjk+XRUcLtyDiSsqkERGQ38WVX/R0QWHOF+fwX8TURSgVXAlbgJapqIXAX8AFx0hPs4KOL3g6r7gns8gGOMMYctFIzJd9TBJIWAiFyC++v9vMi2lCPZqaouBApquWr4kTzu4RC/+1Q0GJsX3BhjDocGY/PD9WC6j64ETgTuU9XvRaQr8LK3YcVOdXNMbVWzMSaBaSgYkzGFA+5BVQuBiQCRVcZZqnq/14HFiqREXgKrf2SMSWQx6uI+mNlHn4hIMxFpBXyNu77gIc8ji5XqloINNhtjEpgGguDxsRTg4LqPmqtqCXAh8JyqDgLO8Das2BFLCsaYJOBOhkmMKan+yLqBi4F3PI4n5moONBtjTMIKBWMyJfVgksLdwD+AFar6pYh0A5Z7G1bsVI8pWP0jY0wi02DI80NxwsENNE8Hpte4vAr4qZdBxVI081pLwRiTwDQUghh0Hx1wDyKSDlwF9AXSq7er6ngP44odG1MwxiQBDQYSY/YR8BLQHhgJfAp0Akq9DCqWbKDZGJMUAkEk5YjWDR+Ug0kKx6rqHUCZqr4AnAPkehtW7EQHmm3xmjEmgWkgkDBJoXoEdoeI9AOaA108iyjGogPNtnjNGJPAYpUUDmbU4qnISuY7cCuZNgXu9DSqGLKBZmNMMkiYpKCqz0T+/BT3qGkNio0pGGOSgQYCSGoCJIXIAXWuwO0yit5eVSd6F1YM2eI1Y0wSSJiWAvAe8G/gWyDsbTixt3vxmiUFY0ziSqSkkK6q/8/zSOJkd/eRDTQbYxKXBgIJczjOl0TkFyKSIyKtqk+eRxYjNtBsjEkGidRSqAIeBG4DNLJNaSCDzjbQbIxJBm5SSPV8PweTFP4f7gK2LV4HExe2eM0YkwQSafHaYqDc60DiZffiNUsKxpjElUjdRyFgoYh8DFRWb2woU1JtoNkYk+g0HHYPx5kgSeHNyKlBsoFmY0yiq+7JiHtSEBEf8HNVbTCH39ybDTQbYxKdVrk9GXEfU1DVEFAuIs09jyReUmyg2RiT2DRQBSRASyGiAvhWRD4Ayqo3NpgxBRFISUGrquIdijHG1Eor3eFcSU/zfF8HkxTejZwaLCctjXBlRbzDMMaYWlUnBSctAZJC5MA6DZqkp6O7LCkYYxJTuLqlkBrHpCAi01T1YhH5lt0rmaupquZ5G1rsOOnp1lIwxiQsrYyMKaTFd0XzpMj5uTW2Ce4xmn/nWURxIOlpaEXlgW9ojDFxoFUJ0H2kqhsi52tEJB+4FLgY+B6Y4XlkMeSkWUvBGJO4ogPN8UwKItIDGANcAmwFpgKiqsPqY8eRNRDzgB9V9dxI5dWpuAfzWQ1crKrb62NfB4wlPd1aCsaYhBXLMYX9rVNYCgwHzlPVk1X1f3BLXtSXScCSGpdvAWarandgduRyTDhpaYQrdsVqd8YYc0hiOaawv6TwU2Aj8LGIPC0iw3HHFI6YiHQCzgGeqbH5fKB6ptMLwAX1sa+DiqdJE2spGGMSVizHFOpMCqr6hqr+F9AL+AT4NdBORJ4QkTOPcL9/Bn7Lnof3bFdjHGMD0La2O4rINSIyT0TmFRUVHWEYLictDa2wMQVjTGIKx3BM4YCls1W1TFX/pqrn4s48WsgRdO2IyLnAZlWdfzj3V9WnVLVAVQvatGlzuGHsGVN6evRFN8aYRLO7+ygBFq/VpKrbgP+NnA7XScAoETkbSAeaicjLwCYRyVHVDSKSA2w+gn0cEifdWgrGmMQV3uUe0sZp0sTzfR3MQXbqlareqqqdVLUL7uymj1T1cuAtYGzkZmOBmbGKSdLSCVtSMMYkqHB5OThOYnQfxdD9wAgRWQ6MiFyOCYm0FFT3XrhtjDHxp+W7cJo0cQt4euyQuo/qm6p+gjuIjapuxZ0CG3NOutsk06qqmGRiY4w5FOHycpyMjJjsK5FaCnHjRMrR2riCMSYRhcvLkQzvxxOgsSaFkg3wr0cJbl1NIBRG0tIBCNtaBWNMAgrv2oWTkQnAtC/Xsu25MTD7bk/21TiTQtlm+OAOrv//nue9bzdED1yhtqrZGJOAwuXl0ZlH975biP/HubDTmwmajTMpNG0HQBspZltZFU66tRSMMYkrvMsdUwiGwpRWVJEZ2gGZ9bNOa2+NMylktEYR2soOtpVVIZGkoFYp1RiTgDTSUijeFaAZ5fg05FlSiOvso7jx+ZHM1nQsL+WrsiqcFu6ofri8PM6BeSMQCLBu3ToqbCDdNHDp6el06tSJlBgc4D6WwuW7cDIy2LErQGspdjdaUqhnTduRU1nCtp1V+Jo2BSBUWhrnoLyxbt06srKy6NKlS0zmORsTD6rK1q1bWbduHV27do13OPWqevbRjvIq2lQnhabWfVS/mralbfWYQlYWAOGdZXEOyhsVFRVkZ2dbQjANmoiQnZ3dIFvE7uyjDLaVBWhL5DAzWTme7KsRJ4V2tNLtbC2rxIm0FMINtKUAWEIwjUJDfJ9rKIRWVOA0yWDrzkraSXVSaO/J/hpxUmhLs9B2tu2s3N19tLPhJoV4Wr16Nf369Tuix9ixYwd/+ctf6imi2Js+fTp9+/bFcRzmzZsX3V5VVcWVV15Jbm4ueXl5fPLJJ9HrXnnlFXJzc+nfvz9nnXUWW7ZsAWDNmjUMHz6c/v37M3ToUNatWxeT5zBv3jwmTpwIwCeffMLnn38evW7cuHG89tprMYmjsQnvcls+TkYGW8uqaCs7UH8TSGvmyf4acVJoR4pWEaooJuzzu+WzG2j3UUOQ7EmhX79+vP7665x66ql7bH/66acB+Pbbb/nggw/4zW9+QzgcJhgMMmnSJD7++GO++eYb+vfvz2OPPQbAjTfeyBVXXME333zDnXfeya233hqT51BQUMCjjz4K7JsUjHfC5e73kpPRhK07q+jo24FktQePWkWNOikAtGEH28urcJo2bdDdR4li1apVDBgwgC+//JK5c+cyZMgQBgwYwJAhQ1i2bBkAixcvZvDgweTn59O/f3+WL1/OLbfcwsqVK8nPz+emm24C4MEHH+S4446jf//+TJ48GXBbJb179+YXv/gFffv25cwzz2TXrl2sXLmSgQMHRuNYvnw5gwYNAmD27NkMGDCA3Nxcxo8fT2Xk2BpdunSJ/jqfN28eQ4cOBeD3v/8948ePZ+jQoXTr1i36RQlwzz330KtXL0aMGMEll1zClClTAOjduzc9e/bc5/UoLCxk+HC35Ffbtm1p0aIF8+bNQ1VRVcrKylBVSkpK6NChwz73GTZsGDNn7r+g8FtvvUV+fj75+fn07NkzOgg7f/58TjvtNAYNGsTIkSPZsGEDAEOHDuXmm29m8ODB9OjRg3/+85+AmwjOPfdcVq9ezZNPPsnDDz9Mfn5+9PrPPvuMIUOG0K1bN2s11KNwWXVSyGBrWSU5vmLPxhOgEc8+2tWkBYvS02gR2si2MncGUrhsZ7zD8txdby+mcH1JvT5mnw7NmHxe3wPebtmyZYwZM4bnnnuO/Px8SkpK+Oyzz/D7/Xz44Yf87ne/Y8aMGTz55JNMmjSJyy67jKqqKkKhEPfffz+LFi1i4cKFAMyaNYvly5czd+5cVJVRo0bx2WefcfTRR7N8+XJeeeUVnn76aS6++GJmzJjB5ZdfTvPmzVm4cCH5+fk899xzjBs3joqKCsaNG8fs2bPp0aMHV1xxBU888QQ33HDDfp/L0qVL+fjjjyktLaVnz55MmDCBr7/+mhkzZrBgwQKCwSADBw6MJp665OXlMXPmTMaMGcPatWuZP38+a9euZfDgwTzxxBPk5uaSmZlJ9+7defzxx6P3mTFjBpMmTeKNN96gtLSUrVu3kp2dXes+Ro0axahRowC4+OKLOe200wgEAvzqV79i5syZtGnThqlTp3Lbbbfx7LPPAhAMBpk7dy7vvfced911Fx9++GH08bp06cJ1111H06ZNufHGGwH461//yoYNG5gzZw5Lly5l1KhR/OxnPzvge8IcWLjE/bz6mjdn6+YqPmy+g+1N2nK6R/trtC2FIr+f8Tnt8GWuYOtOdwZSqLThJ4V4KSoq4vzzz+fll18mPz8fgOLiYi666CL69evHr3/9axYvXgzAiSeeyB/+8AceeOAB1qxZQ5NaDiwya9YsZs2axYABAxg4cCBLly5l+fLlAHTt2jW6j0GDBrF69WoArr76ap577jlCoRBTp07l0ksvZdmyZXTt2pUePXoAMHbsWD777LMDPp9zzjmHtLQ0WrduTdu2bdm0aRNz5szh/PPPp0mTJmRlZXHeeecd8HHGjx9Pp06dKCgo4IYbbmDIkCH4/X4CgQBPPPEECxYsYP369fTv358//vGPAEyZMoVPP/2UAQMG8Omnn9KxY0f8/gP/vvvTn/5EkyZN+OUvf8myZctYtGgRI0aMID8/n3vvvXePsYkLL7xwn9fvQC644AIcx6FPnz5s2rTpoO5jDiwUSQpOs2YU7axgejPhS3/4APc6fI22pdCujfvL1p+ylW1lVXRomtkouo8O5he9F5o3b85RRx3Fv/71L/r2dWO44447GDZsGG+88QarV6+Ods9ceumlHH/88bz77ruMHDmSZ555hm7duu3xeKrKrbfeyrXXXrvH9tWrV5NWo/y5z+dj1y63ptVPf/pT7rrrLk4//XQGDRpEdnY2a9eurTNmv99POOx++Pae5rj3PoLB4GEdj8Pv9/Pwww9HLw8ZMoTu3btHW0THHHMM4P7Cv/9+9xAjHTp04PXXXwdg586dzJgxg+bNm+93P7Nnz2b69OnRhKeq9O3bly+++KLW21c/v+rndjBqviZ2bJL6E9rhrkvwNW9OedkKdrV0aNfEmzUK0IhbCmmZbWgVCiMpJWwrq8TXNKtRdB/FS2pqKm+++SYvvvgi//d//we4LYWOHTsC8Pzzz0dvu2rVKrp168bEiRMZNWoU33zzDVlZWZTWSNojR47k2WefZedO9//sxx9/ZPPm/RcIS09PZ+TIkUyYMIErr7wSgF69erF69WpWrFgBwEsvvcRpp50GuN0k8+e7hxKfMWPGAZ/jySefzNtvv01FRQU7d+7k3XffPeB9ysvLKYv0GX/wwQf4/X769OlDx44dKSwspKioKHpd7969AdiyZUs0Wf3xj39k/Pjx0cfr1avXPvtYs2YN119/PdOmTYu2unr27ElRUVE0KQQCgWhL7WDs/f9hvBMqcZOCZDXDF1gFQNumHTzbX6NNCgDt8VPlL2NTSaV1H8VAZmYm77zzDg8//DAzZ87kt7/9LbfeeisnnXQSoVAoerupU6fSr18/8vPzWbp0KVdccQXZ2dmcdNJJ9OvXj5tuuokzzzyTSy+9lBNPPJHc3Fx+9rOfHdSX1GWXXYaIcOaZZwJuonjuuee46KKLyM3NxXEcrrvuOgAmT57MpEmTOOWUU/D5fAd87OOOO45Ro0aRl5fHhRdeSEFBQfQX/BtvvEGnTp344osvOOeccxg5ciQAmzdvZuDAgfTu3ZsHHniAl156CXBbA5MnT+bUU0+lf//+LFy4kN/97neAO+Dbs2dPevTowaZNm7jtttsAN1nU9gv9+eefZ+vWrYwePZr8/HzOPvtsUlNTee2117j55pvJy8sjPz//kGYTnXfeebzxxht7DDQbb4SK3aSw3Z9GU787GaBt886e7U+SuZlXUFCgNed8H6obXhzCsspS+jR/md+ufI8d01+j51fz6zHCxLBkyZLor8zGbsqUKRQXF3PPPfd48vg7d+6kadOmlJeXc+qpp/LUU0/tMevJS++88w6rVq2KriVorBra+33TH+9nx/TpBN7+iHtevo7C9ov4x+lP0eGoEw/7MUVkvqoW1HZdox1TAGif1oIvQiVsKCnH17Il4fJywlVVOKmp8Q7NeGD06NGsXLmSjz76yLN9XHPNNRQWFlJRUcHYsWNjlhAAzj333Jjty8ROqLgYp3lzNhZX4E/Zik+Vtu36e7a/Rp0UcjLaU16xlm3bNuDr1AqA0LZtOO29WT5u4uuNN97wfB/V4yXG1JdQSQm+5s3ZWFKBppTQLgz+1EzP9teoxxTaNTsagPCu7/C1bAm4ScEYYxJFqLgYX7NmbC6pYFfKLjpK2oHvdAQadVLIadkdgCaylopMt45IcKslBWNM4giXFONr3oyNJRXsSAnRIXX/04+PVKNOCu3b9AEgI2UzW1PcA+2EtltSMMYkjuCWrfhaZbNhRzFbfdAxo62n+2vUSaFN236khRV/6lY2Oe787aB1HxljEoQGAoS2b8ffpg1lxYtREXIi3d5eadRJwfGlcJQ6VKUWsyHoB5+P0Lbt8Q6rwbHS2Q2jdHYsWHnuPQW3bgXAl51NqMJdYNkx0u3tlUadFAA6+5uyI6WS9SWV+Fq1tO6jBJXsSaEhlM6OBSvPvadgZEV7WVYLmqS4yb9jW++mo4IlBTpntGNzCnxftAN/dmuCm4viHVKDZqWz9xSL0tkAEyZMoKCggL59+0Zfq+rnWF0me/DgwdFyH+PGjeO6667jlFNOoUePHrzzzjv7fXwrz+2N6qRQlJKFP3UT6eEw7XO8XfvSqNcpAHRp0Y1g2Qq2b/2GlPbtCUTetA3W32+Bjd/W72O2z4Wf3H/Am1np7H3FonQ2wH333UerVq0IhUIMHz482voAaNasGXPnzuXFF1/khhtuiCaA1atX8+mnn7Jy5UqGDRvGihUrSE9Pr/XxrTy3N4JF7o+S9b4MKlNL6Kx+HF+Kp/u0lkKbXACqyhfjb9+ewMaNcY6oYbLS2bWLVensadOmMXDgQAYMGMDixYspLCyMXnfJJZdEz2tWTb344otxHIfu3bvTrVs3li5desDnY+W561d1S2FFIJVtqVV0TfF2OipYS4GjOxwPCyHF9wNVrXoQLi4mXF6Ok5ER79C8cRC/6L1gpbPr3ofXpbO///57pkyZwpdffknLli2jLaRqNQ92X9fftV3em5Xnrn/BTZvwtWrFii1bKPJDt6xOnu8z5i0FETlKRD4WkSUislhEJkW2txKRD0RkeeS8ZSziyc7uSbNwGCdtM9sy3V1aa6H+Wens2sWidHZJSQmZmZk0b96cTZs28fe//32P66dOnRo9P/HE3UXWpk+fTjgcZuXKlaxatarWMZFqVp7bG4Ef15HSqRPFW+eiInRrVff/QX2JR/dREPiNqvYGTgB+KSJ9gFuA2araHZgduew5cRy6azqlaSWsT3FXNTf4cYU4sdLZ8SmdnZeXx4ABA+jbty/jx4/npJNO2uP6yspKjj/+eB555JE9Wi09e/bktNNO4yc/+QlPPvkk6enprF+/nrPPPnuffVh5bm9UrfuRlI4dqap0u/u6tt//GFW9qJ7lEK8TMBMYASwDciLbcoBlB7rvoEGDtD788dWzteDZvvrwCx9qYc9eum3atHp53ERRWFgY7xASxoMPPqi33367Z49fWlqqqqplZWU6aNAgnT9/vmf72tvbb7+tjzzyyCHdp3PnzlpUVLTP9rFjx+r06dPrK7SYaijv93AwqIX9cvW7e+7Xyx4arvnP9dWKXdvr5bGBeVrH92pcxxREpAswAPgP0E5VNwCo6gYRqXUtt4hcA1wDcPTR9bOyr2d2byoqfmDFrhWMTEkhsGZNvTyuSSxWOtskk+CmTRAIsCUrm8rUrXQLOaSlt/B8v3FLCiLSFJgB3KCqJQcaxKqmqk8BT4F7kJ36iKVXx5Pgx39QsnM+qUcfTeVBzoQwycVKZ++rrlk/Ncd4THxURWZq/ZDSjPVpVQxLi01J/7hMSRWRFNyE8DdVfT2yeZOI5ESuzwH2P2pYj47pejp+VdRZhXY6iipLCsaYOKta9T0Ai4NbKPE55Gb3icl+4zH7SIC/AktU9aEaV70FjI38PRZ3rCEmUtObc0zIoTx9C9tbtSew5ge0xsCnMcbEWuWKFTgZGfwQ+AaAPkeddIB71I94tBROAn4OnC4iCyOns4H7gREishx34DmmE+rzMzvxQ3qQHzJaooEAVT/8EMvdG2PMHipXrMDp2o2gsxy/Kj26jYzJfmM+pqCqc4C6BhCGxzKWmgpyBjN19VqW+DbRB6goLCQtUr/FGGNirXLlCnbmFlCasYze4RTSmsRk6ZaVuag2oIdbt+X71KVIaioViwsPcA+TLMaPH0/btm33Kd/99ddfR9c5nHfeeZSUlAD7L2ddbdSoUb/nqF4AABy+SURBVEdcDvxQ3HnnndHaQH/+858pLy+PXte0adOYxWFiI1hURKhoCz9ktWRNWpiCrNj9QLWkENGu/QBygkpZkzVot2OpOISVlyaxjRs3jvfff3+f7VdffTX3338/3377LaNHj+bBBx8E6i5nXe3111+P+Rfx3XffzRlnnAHsmxRMw7Pr20UALE5ZR0iEE7qcHrN9W1KoJsLgtPasztjFpnZHUVFY2Chqq8TCHXfcwSOPPBK9fNttt/Hoo4+iqtx0003069eP3NzcaLmF6vLJ1f77v/87OkWyS5cuTJ48mYEDB5Kbmxst0lZUVMSIESMYOHAg1157LZ07d46WvT711FNp1arVPnEtW7YsenyDESNGREtZ1FXOGtxaQw899BC33377QT33q6++OlpSuk2bNtx1113AoZX9ht0HmHn00UdZv349w4YNY9iwYXu8pnl5eZxwwgmNolBcQ1ex6FtwHFY1X4lflfxesasG2+gL4tU0rMvpzFz5CsvTd9CutJSq1asb3LjCA3MfYOm2A1e7PBS9WvXi5sE313n9VVddxYUXXsikSZMIh8O8+uqrzJ07l9dff52FCxfy9ddfs2XLFo477rh9DkJTm9atW/PVV1/xl7/8hSlTpvDMM89EC93deuutvP/++zz11FMHfJx+/frx1ltvcf755zN9+vRocbz9lbO+4447+M1vfkPGQRZMfOaZZwC3NtDIkSMZN27cYZX9rjZx4kQeeughPv74Y1q3bg1AWVkZJ5xwAvfddx+//e1vefrppw86aZnEtGvhQoJHd2VD8zXkhlLJyIrNGgWwlsIeTug/Dr8qi1utBqD83/+Ob0ANRJcuXcjOzmbBggXRctfZ2dnMmTOHSy65BJ/PR7t27TjttNP48ssvD/h4tZVcnjNnDmPGjAHgrLPOomXLAw/KPfvsszz++OMMGjSI0tJSUlNTgbrLWS9cuJAVK1YwevToQ3r+FRUVXHTRRTz22GN07tz5sMp+709qamq0ZXUoZahNYgpXVlI+/yvWdGjL2jQY1mZATPdvLYUaMpt1IC+URmFOEaG27Sj7/HNaRmrNNxT7+0Xvpauvvprnn3+ejRs3Rqt61tU9V7NkNdRdtrpmyeXD6err1asXs2bNAuC7776LVjWtq5z1p59+yvz58+nSpQvBYJDNmzczdOjQWgeia7ruuuu48MILo2MCehhlv/cnJSUlWtb6UMpQm8S0a8ECtLKSb1qsB+CMvHEx3b+1FPZydseT2ZDmsLpjNmVf/Bu1D1i9GD16NO+//z5ffvlltELoqaeeytSpUwmFQhQVFfHZZ58xePBgOnfuTGFhIZWVlRQXFzN79uwDPv7JJ5/MtGnTAPcAPNu3bz/gfapLbYfDYe69995oddS6yllPmDCB9evXs3r1aubMmUOPHj2iCeGNN96o9VjJjz/+OKWlpdxyy+6iv4dT9rsmKyndsJX963Pw+VjYZT3dgnBU51Niun9rKexl5PG/5o8zZ/Ntmw0cs2AnuxYuJKOgIN5hJb3U1FSGDRtGixYtomWoR48ezRdffEFeXh4iwp/+9Cfat3f7Ti+++GL69+9P9+7dGTDgwM3nyZMnc8kllzB16lROO+00cnJyyMrKAtwjin3yySds2bKFTp06cdddd3HVVVfxyiuvRA9xeeGFF0aPsbB582ZGjhyJ4zh07NgxWs56f1auXEmzZs322T5lyhRSUlKiXULXXXcd1113HUuWLIkeu6Bp06a8/PLLB1WeG9yiez/5yU/Iycnh448/Pqj7mOSx85//ZOtRHVnSfD2/zMqPfQB1lU9NhlN9lc7e27V/PV6H/6W3LurXXzfcdbcn+4ilRCglHAqFNC8vT7/77jtPHr+iokIDgYCqqn7++eeal5fnyX7qctlll+nmzZtjuk9Tu0R4vx+uihUrtLBnL33sqiHa/7m+umn9Ak/2w35KZ1v3US0u73sJmzJ8fN+1OSXvv48GAvEOKakVFhZy7LHHMnz4cLp37+7JPn744QeOO+448vLymDhxYnStQay8/PLLtGnTJqb7NA1Pybvvoo7D33O3c1wonbY5sW8pWPdRLYYU/JKjv32af/baxjFvhdj5zzlknT7swHc0terTpw+rVq3ydB/du3dnwYIFnu7DGC9pKETxW2/zY8cs1rQo45ZeYw98Jw9YS6EWjs/PpTnD+EcvKM9KZ9tLL8Y7JGNMA1f60UcE1q3jg/xSjgnAScf9Mi5xWFKow8Uj7qe9wsd5VZR/8W8qliyJd0jGmAZKVdn2/AuUNkvlH32Ea48dgzjx+Xq2pFCHlNRMrjrqXF4rEKpSfRQ98mi8QzLGNFBlc/7Frvnzeee4AH1Cfs46dd/pzbFiSWE/fjr8Hrr5/bx1Qpidn3xC2b//E++QjDENjAaDbHrwQbY3c3h7oHDTcb+NWysBLCnsl+NLYfJpD/JugbC9mfDj7bcTjixqMkfu6quvprCwfkqUx7t89Pr16/nZz9yiZQsXLuS9996LXvf73/+eKVOmxCs0k+C2Pv00Vd99x1/PgPNTuzAw79K4xmNJ4QB6dh/BFdnH8+fzHILr1rHxgT/FO6QG45lnnqFPn9gcd9ZrHTp04LXXXgP2TQrG1GXX4sVsfuxx/tMLiro53HrR3+IdkiWFgzFh9DO069Cct48XiqdNY9urr8Y7pKRSVlbGOeecQ15eHv369YuWyB46dGi0JHXTpk25+eabGTRoEGeccQZz585l6NChdOvWjbfeeguA559/nvPPP5+zzjqLnj17RstQ7622stR1mTdvXrS0dW5ubrSG0MqVKznrrLMYNGgQp5xySrRE97hx45g4cSJDhgyhW7du0USwevVq+vXrR1VVFXfeeSdTp04lPz8/+lwLCwujz+fRR218ykBg82ZWXXMNxU1CvDJCeGjY46Slt4h3WLZO4WCI4/DgJW9zdfB0jtpSSd5dd+NvlU2zM0fEO7RDtvEPf6BySf2Wzk7r3Yv2v/tdnde///77dOjQIVpwrri4eJ/blJWVMXToUB544AFGjx7N7bffzgcffEBhYSFjx45l1Cj3yHhz585l0aJFZGRkcNxxx3HOOedQUKMMSV1lqesqyV1QUMDChQsBuOmmmzjrrLMAt5TEk08+Sffu3fnPf/7D9ddfz0cffQTAhg0bmDNnDkuXLmXUqFHRbiNwy3ncfffdzJs3j8ceewxwu4+WLl3Kxx9/TGlpKT179mTChAmkpKQc9GtsGpbg9u0sGzeWYMk2HrzM4ZaB/82xxxy4bHwsWFI4SOkZrXjykveZEBxJ5owqmDSRtvfcS5uf/TTeoSW83NxcbrzxRm6++WbOPfdcTjll3wJfqamp0S/k3Nxc0tLSSElJITc3d49S0CNGjCA7Oxtw6xXNmTNnn6RQXZYa3IPiLF++/IDHaZg2bRpfffUVs2bNYufOnXz++edcdNFF0esrKyujf19wwQU4jkOfPn0O+oA255xzDmlpaaSlpdG2bVs2bdpEp06dDuq+pmEJbNzIokv/C9+mzTxyocPVJ/6CYSdOiHdYUZYUDkHTZu15evwn3Cg/oXJmCc7tt1P05Rf0vucPSKQWf6Lb3y96r/To0YP58+fz3nvvceutt3LmmWdy55137nGbmuWfHceJlpB2HGePUtDVt6nrstZRlnp/Fi9ezOTJk/nss8/w+XyEw2FatGgRbUHsrWZ5az3Ikt17l8S28taNU/G/PmfVxOvQQIDHfypcee7NnDr4iniHtQcbUzhE6U1a8ui1c1jz8zzeO06Qme/ynxEnU7Kg9i8Q487MycjI4PLLL+fGG2/kq6++OuzH+uCDD9i2bRu7du3izTff5KSTTtrj+v2VpR4+fDg//vjjHrcvLi5mzJgxvPjii9HaRc2aNaNr165Mnz4dcL/4v/7664OO0Upbm72Fdu5k/qRr+fGqq9iSHuCvl6Zw2zUvJVxCAGspHBbH5+emy19hTp9pvPTq3ZzzUSlrL72EooHdGHzPY2R0a1iH8DxS3377LTfddBOO45CSksITTzxx2I918skn8/Of/5wVK1Zw6aWX7tF1BHDmmWfWWpa6devWrFixYp9jNb/55pusWbOGX/ziF9FtCxcu5G9/+xsTJkzg3nvvJRAIMGbMGPLy8g4qxmHDhnH//feTn59f6zEWTOOhVVV8+8TDVL74IhllYT4cKATPzOeRy57Dn5J24AeIAznY5m8iKigo0OrZK/ESDgV46tWJVL3/CactAF8YNvRtR48Jt3L0sDP36d6IhyVLltC7d+94h3HEnn/++T0GcA/FokWLePbZZ3nooYc8iMwkkkR4v1du2swXf55Mkw8/o1lpmO86wNfDWnLZ2Mc55ujYHl6zNiIyX1VrPVCMtRSOkONL4brLnqDiwhJeeOVGdPYcjlu0ifLrb2BOtkPZoL7kXXkjOQMGxzvURq1fv36WEIynAsXFzH3+YUpnv0fHFaW0C8PSo+CfI7M58/J7OL9XclRatpZCPdNwmH989ATfz3yBjoWldI90YW9pKezo1prMwSdT8F+/pFn7jjGLKRF+ORkTK7F6v2sgwOJZM1g1axqpS1fQYV2AlBDsyIQlPXwEh+Qz+rL7aNOqs+exHKr9tRQsKXiofFcx77/zZ7Z9/D4tv99Bt3WQHoAwsLWlUNyuCcHOnWiZfwJ5Z46hZUdvxiIsKZjGxIv3e7CiksLPZrL237OoWr6MjA3babcpRFrk+FvrW8G6rqkE83pz8k9/Tc9jjq/X/dc3SwoJYul3X7DgvWcJfPsNGRtLydmstKoxSaW0Cexo6VDWKp1A25akduxC62Ny6Zo/lJzu/XAO8hi+e1uyZAm9evVKiPENY7ykqixduvSwkkJJ0UaWf/UJG7+bR8WaFbBxM+lbd5K1I0DLYvCH3dtV+WF9a9jeLo3KYzrS4eSfMHT4laSlZdbzs/GOJYUEVVlVzoIv3uKHz/9B5ZoVpGwpJmt7gFY7oMVedfeCDpQ0hdIsoSLTTzAzjXDTTJxmLfG3ak1GmxyyOnalbeeedOjclyZZzaP3/f7778nKyiI7O9sSg2mwVJWtW7dSWlpK165dCQUCFK1dzqbV37F9/UpKN62hcstGQiXb0Z07SS3dRXppgMyyMFk7oUnVno9X6YeiFlDSwqGyVSbaIYfm/QrIH3YRnTr2is+TrCeWFJJMOBRi9fIFLJ//IdtXFxLY/CO+7aWklVSQURqkSbnStBwyK+t+jMoUqEiFyjQItGpO08uuJSXnKHAcVAARQEDETRSOgDhI9ckRxPEhjoM4Phyn+jpLKqb+aVgJawgNhQiHw6iG0HAY1TAaDoMqqu757hNI5BwUCYcJrl9LxfP/S/qWEjJ31b0QKySRH1lNoSzTR2VWCqFmTfG3bktmh26073cCvQaeTlZWyxi+CrGTVElBRM4CHgF8wDOqen9dt22oSeFgqCpFRT+w9rv5bP9xFWUb11K1dTOh4u1oWRmyqwKnMoC/KkhKZZjUSiWtCtIr3VPaYSyoDQkE/BByIOiDUOQUrv7bkchlIeQTwnud1HHAiZz7HHBqnnyIzxc9d2/jJiX1VV/nID6/u93n/i0+P47Ph4iD4/Mjjg+fz+cmMp8fcfyRpOaei8/v3tbxIY64j+n4EBG3hn3kXBwBiSRFEXd/kYTpJlJwVNzkSeQ+qruTrQIiqArut5egYffLTMNhlDBa/beGUFXCVH8Bul+GYbTGF2MI1P3ilGCYkIYgHIJwePcXaSiEhqu/WJVwyL2NezmEBgOEQ0E0HCQcCqPhoHv76Ln7RVz9mBJyH59QEAmFIRSOnIeQsEIojBNyv4wlHEZCCuEwElacsEbOQcKKhMEXUvxB8AchJRQ5j5wOZxVtRQpUpkJVinsKpArBVMc9pacQzEhDM5ogzZrhb9GSjDadyO54LB2OyaNjl774/I138mXSTEkVER/wODACWAd8KSJvqWr9FN1vQESEtm0707bt4c1sqKrYxdaidewo+pHiLevZVVxExc6tVO4sIVxeRqiinHBFBVRWoJVVSFUVBIIQDLpfFqEQTtD9QnAiXw5OyP3g+wJKamUYX4g9Tv4QOAqi7noOJ7y7n9Ykj6Dj/jAIOdU/BiDs7D7XyN/qCGEfqLjnlWkOu3wOIb8Q9juE/A7q96Ep1Sc/kpqOLzUdSW+CPy2DlCaZpGQ0Jy2zOekt25DRog3N2nYku10nspq2jOvBaBqqhEoKwGBghaquAhCRV4HzAUsK9Sw1vQk5R3Un56jucYtBVQmFwlRUlhOoKKWiooyqijICFWWEgkFCVRWEApUEAwFCgUpCgQDhYIBQoIpQIIAGA+62cJBwKEQ4HERDIUKRX77hcOSXcjgU7Yao7pYQDaMKhBWi3RJEuik00i2xZ1cFGkbDgAC429yGgBJtb8vulrdG/hUit3Nv4N69utuuugsvep1EWxsi1dft/lvFibRonGiXHzVbPNWtHJ97O8fxgd+P4/jx+VJw/H4cXwo+fwp+f6r7d0pq5HI6Pn8qfn8qvrQ0UlKakNIkA39aOilpmaSkZ5CenkFKit+6ERuwREsKHYG1NS6vAxJ7bpc5bCKC3++jqT8LMrPiHY4xhsRLCrX9/Nhj0ENErgGuiVzcKSLLPI9qT62BLTHe5+GyWL1hsXrDYvVGbbHW2e+caElhHXBUjcudgPU1b6CqTwFPxTKomkRkXl0DNInGYvWGxeoNi9Ubhxproo3SfAl0F5GuIpIKjAHeinNMxhjTaCRUS0FVgyLy38A/cKekPquqi+McljHGNBoJlRQAVPU94L14x7Efceu6OgwWqzcsVm9YrN44pFgTbvGaMcaY+Em0MYWEIzYh25h6Z5+rxGVJoRYikisiNwNoAjelRKSziBwd7zgOhojkJMsXgYi0i3cMB0tEeibReyApPlfQuD9blhRqENcU4P8Av4ikxDum2ohIExF5GPg78IKITIhsT7j/TxFJE5EngE+Bp0TkwnjHVBcRaVr9uorI/yZyrAAi0hJ3tf/VItI63vHUJVk+V2CfLbCksLc2QA4wSFXvU9VAvAOqw0Sgg6r2AX4P3ACgqolYSWgUkKOqPYB3gLtFpEecY9qHiHQEXsL9TJyN+0H7U1yDOrCOwFIgE8iPcyz7kyyfK7DPliUFEalZX6E50F1Vq0RkpIjcKCIj4xVbTSLSNHLu4P6/fRO5qgPwrogkVIH3SHFDgBBQBKCqM4H3gWtFpEW8YqtDBW5V3kmquhGYBiwUkf5xjmsPkV/d1V0FxcBruKv+h4lIdvwi21eNX9fNSNDPVTUR8UXiFRL4syUiGTUuKh58thptUhCRY0VkGvCciJwjIq2AncC/RORu4Le4XxR/FpGx1V/KcYzzBRE5B3f9xlKgm4j8E3gAKAU+FJER8ey3j/RvPykiTVQ1FNmcCmwTkQ6Ryw8CA4G+kfvEJd6asQKo6lbgkxo3OQroBsS6jMo+9npdtUZ//ECgCXA70Ba4REQuEJH0BIm1+te1A3yWSJ+rWmINReJdDxydgJ+t7iLyIm430ajI+7YK2FHvny2NVIVsTCfcN+k7wB24VVj/AkyJXPc/uF8OeZHLP8P9NZaVAHE+CTwYua457q/Z9pHLvwTei+NrejIwF/cQ1LfV2J4TeQ5nAamRbb8HXkukWIlMz65xm57A6/GK8QCxOpHzY4BJkb8/wm2V3ZJIsUa2twQeBT5OhM/VAWJ1EvCz9XPcsaMJwHjgr8BPgVZefLYaa0shB9gB/EHdZtd9wCkicirwDBAABgGo6mtANhCPMp57x3kPcLKIjMJt1WwCqpuTzwAZcexC2Ir7hu0BXCkiXQBUdQPuCvUL2d3v/SqwPY4DjnvH2lkjn6gaBgArAUTkF3HsRtrnddXdv8BPxB1kXgRsBF4B1u7VxRBLdb0HtuN+eVWQGJ8rqCNW3C6ZpsCPJM5naxNws6o+oarP4ha3a6+q24APqOfPVqNMCqr6I1CAezCf6stPAHeo6te4KwDPE5FbI83IRcC2BIrzV+p2zxwFXCMi43C/eL/E7WeOOVVdgnssjBW4b9S7alz9F9yS6LeIyG9w37irNE4DjrXEejfsM8NkOJAtIjOAS3G/0GLuAK/rq7j/5+NV9VLgTdz3RFy6OQ4Q6wfAi8C58f5cQd2xRn4cbAS6A79IkM/WLGCWiFRXoKjAHe8AeIz6/mzFq0kUo2ZXDjBkr22+yPk4YE6N7S2AqdW3x+2X+w0wJgHjnB6Jryduk/LdWMRZV6w1rqteIZ8FrACG17guHTgN91CrlydyrJHtfwcWAz9L9Fhr3M5J9FiB3Fh9ro4kVqA/cH2ifLb2ut3fgAtrXE4HTq2vz5bnTzSep8iH+u/AwJpvgsjfPtx+2BtqbHsB6JcEcb4YjzgPFGt1vJHzG4B3In9fAmQkUawCDE2SWNOT6HVNpvdAkwSM1cHt0noDaBd5n44E0uozjgbZfSQiTmRWw1rc6WWniUimqmp1F4G63S83AZMiMzYuB/rhDjwlepx9YxnnQcRas7siHIn7z8BJIlIMnB55jJh0axxhrGfgfsg+SYJYTyeGXcCN6D1weiTUhIpV3bGk5pHTOcB/gFPqO9YGkxSkxtxcVQ2r6k7cJmER0IXdA1zhyO19qjof9xfCANyjuU1SVU+PB50scR5irFrjdioizUXkQdw3+dmq+gtVLa95uwSO9SpV9XT8oL5f12SKNUneAwkZa0QBMBR3xtFvVPV2Va2o11jrs9kRjxPuDIbHgFnAr4Deke3HAs9H/v5/uHOObwA6aS1NM4vzyGOtcX8H6G+xWqwWa73Hmglc62WMSd1SEJETcPvgtgD34h6+83oAdWcVlIi7srYnMAk4XlXXRa737FdAssZ5pLFWU/eXzzd4zGK1WBtTrCLiqGqZqv6vl3Em3EF2DtEO4CFVfRWiBcKGR+bo9sGd6fAN7nzzKYCKSHdVXW5xWqwWq8WaTLFqjOovJV1SEBGp/vWsqktFZG2NbQHgWHXn6H4tInOAf6rqLBHpinvM55jMN0+WOC1Wi9VitVj3iDnGvROHRdwVmqOAmaq6az+3ux44WlVviVlwe+4/KeKMxGCxesBi9YbFGjsJP6Yg7orCRbjLuDNrbHdq+bsrbj0TRORCEWlvce7LYvWGxeoNizW2Err7SNxaI2cD5+peUzB195TN/jUGifoCW0Tk77jL5+dYnBarxWqxWqyHQGMwDetQToC/xt/tgbcif3cDrmL3ar92uCt75+DWATkaKMGtcHqBxWmxWqwWq8V6GM8p3gHUfHFxR9wfAc6IbOuBW/nxPODfuHN3F+Eu4DoetzBczce40uK0WC1Wi9ViPYLnFu8AIi+O4FbSfBm4DPgQmBC5bhowGxgRuXwy8D2QUuP+qRanxWqxWqwWaz08v3gHEHmRmgGfEzngBm6Rp8dxl3L3BJYDZ9W4/TtArsVpsVqsFqvFWr+nhJh9pKolwGrcMtEA/wLm4R5tbCXwLDBC3IOd/BV31tRKi7NuFqs3LFZvWKyJIyGSQsQbQL6I5KhbHOpr3MUdxwAP4zbV+gMrVfVs9bgQWAOIEyxWr1is3rBYE0G8myo1mlg5wJ+AW2tsmwOcVOOyz+K0WON9slgt1mSK9VBPCdNSUPdYvm8CPxGRi8Q9ZmoFbvatvk0oPtHtlixxgsXqFYvVGxZrYki4Mhci8hPgImAI8JiqPhbnkGqVLHGCxeoVi9UbFmt8JVxSABC3aqCqajDesexPssQJFqtXLFZvWKzxk5BJwRhjTHwkzJiCMcaY+LOkYIwxJsqSgjHGmChLCsYYY6IsKRhjjImypGCMMSbKkoIxxpgoSwrGGGOi/n9H6yBuNsQgZgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "pvlib.atmosphere.get_relative_airmass(solpos['zenith']).plot(label='kastenyoung1989, zenith')\n",
    "pvlib.atmosphere.get_relative_airmass(solpos['apparent_zenith']).plot(label='kastenyoung1989, app. zenith')\n",
    "pvlib.atmosphere.get_relative_airmass(solpos['zenith'], model='young1994').plot(label='young1994, zenith')\n",
    "pvlib.atmosphere.get_relative_airmass(solpos['zenith'], model='simple').plot(label='simple, zenith')\n",
    "plt.legend()\n",
    "plt.ylabel('Airmass')\n",
    "plt.ylim(0,100);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEKCAYAAAAFJbKyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3deZgcVb3/8fe3pyczmZmQjWwQJAECgQDZBhUQDAZCrkAgIoJeMAiK6zXqLwh4VVB/8eKFex/1QdBcEPJTlFUgIGBC2IIoYbJAyMJN2MKQkMlCtklm//7+qJqZzjCT6c50d/X0fF7P0091ndq+VZD+zjlV55S5OyIiIsmKRR2AiIh0L0ocIiKSEiUOERFJiRKHiIikRIlDRERSosQhIiIpyVjiMLPfm1mVmb2WUDbAzBaY2dpw2j9h2XVmts7MXjezszMVl4iIdE0maxx3AVPblF0LLHT3UcDCcB4zOw64BBgTbnOrmRVkMDYRETlAGUsc7v48sK1N8fnA3PD7XOCChPJ73L3W3d8C1gEfzVRsIiJy4OJZPt4Qd98I4O4bzWxwWH4o8M+E9SrDsg8xs6uAqwBKS0snjh49OoPhioh0TV1DE69v2kVhQYzRQ/tEHQ4AS5Ys2eLugw50+2wnjo5YO2XtjoXi7nOAOQDl5eVeUVGRybhERLrkjc27mfxfzzHy4FKemTUp6nAAMLN3urJ9tp+q2mRmwwDCaVVYXgkclrDecGBDlmMTEUm7hsbgb+B4rL2/j7unbCeOecCM8PsM4JGE8kvMrMjMRgKjgMVZjk1EJO3qG5sAiBfkT++HjDVVmdmfgUnAwWZWCVwP3AjcZ2ZXAuuBiwDcfaWZ3QesAhqAb7p7Y6ZiExHJloamoMbRqyB/ahwZSxzu/vkOFk3uYP3ZwOyuHre+vp7Kykpqamq6uqu8U1xczPDhwyksLIw6FJEeo0E1jtxXWVlJnz59GDFiBGb5k+G7yt3ZunUrlZWVjBw5MupwRHqMet3jyH01NTUMHDhQSaMNM2PgwIGqiYlkWfM9jsI8qnHkz5kkUNJon66LSPY1NDU3VeXPv7+8TBwiIrmitakqf35u8+dMctyXv/xlVq1alZZ9lZWVpWU/IpJ5zf04CvOoxpF3N8dz1e233x51CCISgdamqvz5Oz1/ziSHVFdXc8455zB27FiOP/547r33XiZNmkTz8ChlZWVcc801TJw4kTPPPJPFixczadIkjjjiCObNmwfAXXfdxfnnn8/UqVM55phj+MlPftLusW666SZOOukkTjzxRK6//vqsnaOIJKdeNY7u5SePrmTVhp1p3edxhxzE9eeN2e86Tz75JIcccgh//etfAdixYwe33XZby/Lq6momTZrEL37xC6ZPn84Pf/hDFixYwKpVq5gxYwbTpk0DYPHixbz22muUlJRw0kkncc4551BeXt6yn/nz57N27VoWL16MuzNt2jSef/55Tj/99LSes4gcuOZ+HIW6xyH7c8IJJ/DUU09xzTXXsGjRIvr27bvP8l69ejF16tSWdT/5yU9SWFjICSecwNtvv92y3llnncXAgQPp3bs3n/nMZ3jhhRf22c/8+fOZP38+48ePZ8KECaxZs4a1a9dm/PxEJHn1Yc/xfHqqKq9rHJ3VDDLl6KOPZsmSJTz++ONcd911TJkyZZ/lhYWFLY/GxmIxioqKWr43NDS0rNf28dm28+7Oddddx1e/+tVMnIaIpEF9g/pxSBI2bNhASUkJl156KbNmzWLp0qUHtJ8FCxawbds29u7dy8MPP8ypp566z/Kzzz6b3//+9+zevRuA9957j6qqqvZ2JSIRabk5nkc9x/O6xhGVFStWcPXVVxOLxSgsLOS2225j1qxZKe/nE5/4BJdddhnr1q3jC1/4wj73NwCmTJnC6tWrOfnkk4Hgpvsf//hHBg8e3N7uRCQCLf048qjGocSRAWeffTZnn332PmXPPvtsy/fmGgLADTfcsM96icsGDx7MLbfc8qH9J64zc+ZMZs6c2cWIRSRT8rEfR/6kQBGRHNTaVJU/P7eqceSoyy+/nMsvvzzqMESki1r6ccRV4xARkSSoH4eIiKSkIQ/7cShxiIhkUF0evgEwf85ERCQHtTZVqcYhKbjhhhu4+eabO1y+efNmPvaxjzF+/HgWLVqU8v7vuusuvvWtbwHw8MMPp234dhHpuoY87MeRP2fSjS1cuJDRo0ezbNkyTjvttC7tS4lDJLc0j1WlfhzSqdmzZ3PMMcdw5pln8vrrrwPwxhtvMHXqVCZOnMhpp53GmjVrWL58Od///vd5/PHHGTduHHv37uXrX/865eXljBkzZp+h0keMGMGWLVsAqKioYNKkSfsc88UXX2TevHlcffXVjBs3jjfeeCNr5ysi7WvIw3eO53c/jieuhfdXpHefQ0+Af7lxv6ssWbKEe+65h2XLltHQ0MCECROYOHEiV111Fb/97W8ZNWoUL730Et/4xjd4+umn+elPf0pFRUVLL/HZs2czYMAAGhsbmTx5Mq+++ionnnhip6GdcsopTJs2jXPPPZfPfvazaTldEema1lfH5k+NI78TR0QWLVrE9OnTKSkpAWDatGnU1NTw4osvctFFF7WsV1tb2+729913H3PmzKGhoYGNGzeyatWqpBKHiOSe5p7jqnF0F53UDDKp7RDoTU1N9OvXj+XLl+93u7feeoubb76Zl19+mf79+3P55ZdTU1MDQDwepyn8n7C5TERyW+vN8fypceRPCswhp59+Og899BB79+5l165dPProo5SUlDBy5Ejuv/9+IHiXxiuvvPKhbXfu3ElpaSl9+/Zl06ZNPPHEEy3LRowYwZIlSwB48MEH2z12nz592LVrVwbOSkQOREs/DvUcl/2ZMGECF198MePGjePCCy9seVLq7rvv5o477mDs2LGMGTOGRx555EPbjh07lvHjxzNmzBiuuOKKfd7Bcf311zNz5kxOO+00CgoK2j32JZdcwk033cT48eN1c1wkB7TeHM+fGoe5e9QxHLDy8nKvqKjYp2z16tUce+yxEUWU+3R9RLJr2i0v8GrlDh7+5qmMO6xf1OEAYGZL3L288zXbpxqHiEgG1et9HCIikop87MeRP2eSoDs3v2WSrotI9rWMjptH/TjyLnEUFxezdetW/Ui24e5s3bqV4uLiqEMR6VHqGvKvxpF3/TiGDx9OZWUlmzdvjjqUnFNcXMzw4cOjDkOkR2l5dWwe3ePIu8RRWFjIyJEjow5DRARI6ACofhwiIpKM+jzsxxFJ4jCz75rZSjN7zcz+bGbFZjbAzBaY2dpw2j+K2ERE0qm553iveP78nZ71MzGzQ4FvA+XufjxQAFwCXAssdPdRwMJwXkSkW2u+Od4rj26OR3UmcaC3mcWBEmADcD4wN1w+F7ggothERNImfBpXbwDsCnd/D7gZWA9sBHa4+3xgiLtvDNfZCAxub3szu8rMKsysQk9OiYhkXxRNVf0JahcjgUOAUjO7NNnt3X2Ou5e7e/mgQYMyFaaIiHQgirrTmcBb7r7Z3euBvwCnAJvMbBhAOK2KIDYREelEFIljPfBxMyux4G1Hk4HVwDxgRrjODODDY46LiEjkst4B0N1fMrMHgKVAA7AMmAOUAfeZ2ZUEyeWijvciIiJRiaTnuLtfD1zfpriWoPYhIiI5LH+eDxMRkaxQ4hARkZQocYiISEqUOEREJCVKHCIikhIlDhERSYkSh4iIpESJQ0REUqLEISIiKVHiEBGRlChxiIhISpQ4REQkJUocIiKSEiUOERFJiRKHiIikRIlDRERSosQhIiIpUeIQEZGUKHGIiEhKlDhERCQlShwiIpISJQ4REUmJEoeISIY0NXnUIWSEEoeISIbUNTYBUBTPr5/a/DobEZEcUlPfCEBxYUHEkaSXEoeISIbUNgQ1juLC/Pqpza+zERHJIc01jqK4ahwiIpKEmnrVOEREJAV7wxpHb93jEBGRZOypbQCgdy8lDhERSUJ1XVDjKO0VjziS9FLiEBHJkOqwxlFWrMQhIiJJ2B0mjtIiJQ4REUlCc42jjxJH15lZPzN7wMzWmNlqMzvZzAaY2QIzWxtO+0cRm4hIuqjGkV6/Ap5099HAWGA1cC2w0N1HAQvDeRGRbmtXjRJHWpjZQcDpwB0A7l7n7tuB84G54WpzgQuyHZuISDqpqSp9jgA2A3ea2TIzu93MSoEh7r4RIJwObm9jM7vKzCrMrGLz5s3Zi1pEJEVqqkqfODABuM3dxwPVpNAs5e5z3L3c3csHDRqUqRhFRLpstx7HTZtKoNLdXwrnHyBIJJvMbBhAOK2KIDYRkbRp6cdRpJ7jXeLu7wPvmtkxYdFkYBUwD5gRls0AHsl2bCIi6ZSvTVVRnc2/AXebWS/gTeBLBEnsPjO7ElgPXBRRbCIiaVFdGww5UtbTEoeZHUnQtFRrZpOAE4H/Fz4JdUDcfTlQ3s6iyQe6TxGRXLOrph7Iv8SRTFPVg0CjmR1F8AjtSOBPGY1KRCQPtAxy2AMTR5O7NwDTgV+6+3eBYZkNS0Sk+2tscgAKC/JrdKdkzqbezD5PcMP6sbCsMHMhiYhILksmcXwJOBmY7e5vmdlI4I+ZDUtERHJVpw1v7r4K+DZAOPBgH3e/MdOBiYhIbuq0xmFmz5rZQWY2AHiFYKiQ/858aCIikouSaarq6+47gc8Ad7r7RODMzIYlIiK5KpnEEQ+HAPkcrTfHRUSkh0omcfwU+Buwzt1fNrMjgLWZDUtERHJVMjfH7wfuT5h/E7gwk0GJiEjuSmbIkWLgSmAMUNxc7u5XZDAuERHJUck0Vf0BGAqcDTwHDAd2ZTIoERHJXckkjqPc/UdAtbvPBc4BTshsWCIikquSGnIknG43s+OBvsCIjEUkIiI5LZkhG+eEPcZ/RPCypTLgxxmNSkREclYyT1XdHn59Djgis+GIiEiuS+apqn7AFwmap1rWd/dvZy4sERHJVck0VT0O/BNYATRlNhwRkfywN3yJU1E8v97FAckljmJ3/17GIxERySPb9tQB0L+kV8SRpF9S/TjM7CtmNszMBjR/Mh6ZiEg3tr05cZTmX+JIpsZRB9wE/DvgYZmjG+UiIh3avifoydC/JP9emJpM4vgeQSfALZkORkQkX3wQ1jj65WHiSKapaiWwJ9OBiIjkkw9aahw9s6mqEVhuZs8Atc2FehxXRKRjH1Tn783xZBLHw+FHRESS1HyPIx+bqvabOMysALjM3fWqWBGRFGzvqY/junsjsMfM+mYpHhGRvLAtj2+OJ9NUVQOsMLMFQHVzoe5xiIh0bMvu4JbwwWVFEUeSfskkjr+GHxERSVLVziBxDDmouJM1u59kRsedm41ARETySdWu5hpH/t3j6DBxmNl97v45M1tBa4/xZu7uYzMbmohI9xcv6FmDHM4Mp+cmlBnBO8d/kLGIREQkp3WYONx9Yzh9x8zGAV8APge8BTyYnfBERCTX7K+p6mjgEuDzwFbgXsDc/YwsxSYiIjlof01Va4BFwHnuvg7AzL6blahERCRn7e+uzYXA+8AzZvY/ZjaZ4B5HWphZgZktM7PHwvkBZrbAzNaG0/7pOpaIiKRPh4nD3R9y94uB0cCzwHeBIWZ2m5lNScOxZwKrE+avBRa6+yhgYTgvIiI5ptPnxNy92t3vdvdzCZ6oWk4Xf9TNbDhwDnB7QvH5QHOfkbnABV05hoiIZEZKDxi7+zZ3/527f6qLx/0l8H2gKaFsSMKTXBuBwe1taGZXmVmFmVVs3ry5i2GIiEiqst4zxczOBarcfcmBbO/uc9y93N3LBw0alOboRES6rr6xqfOVurFkxqpKt1OBaWb2aaAYOMjM/ghsMrNh7r7RzIYBVRHEJiLSZe/vqAHgkL75N04VRFDjcPfr3H24u48g6CfytLtfCswDZoSrzQAeyXZsIiLp8O624G3bwweURBxJZuTSICo3AmeZ2VrgrHBeRKTbefeDIHF8JE8TRxRNVS3c/VmCR31x963A5CjjERFJh3e37QXgsP75mThyqcYhIpIXmmschw3oHXEkmaHEISKSZuu3NSeOdmocy+6GFQ9kOaL0irSpSkQkHzU3VbV7j+ORbwTTEz6bxYjSSzUOEZE0a37f+KA8fN84KHGIiGRMLNZmXNja3dEEkmZKHCIi2fL2omB6xKQoo+gyJQ4RkWxZ91QwPbJ79zxQ4hARyZa184PpqLOijaOLlDhERLJl+/pgOmh0tHF0kRKHiEgaNTV55ytZ2l6mGgklDhGRNNqwI+jDMfSg/BwZF5Q4RETSal1V8MjtUYPL9l2wZ1swte7/s9v9z0BEJId0mDj+98lgOvrcLEeUfkocIiJptGrDTgBGD+2z74LVjwXTY8/LckTpp8QhIpJGr1RuB2DsYf32XfD6X4PpqClZjij9lDhERNLojc3VAIxq21TVrHe/9su7ESUOEZEMiBfk789r/p6ZiEiuqNkZdQRppcQhIpJprz8RTI85J9o40kSJQ0QkTTrsNb7ivmDajV/elEiJQ0QkTSo/6KDXePOouMd8OssRZYYSh4hImlS8E/QOnziif/srFObHMCRKHCIiabLknQ8AKD88IXF4EoMedjNKHCIiadKaOAa0Fq7/RzAdPCaCiDJDiUNEJE3WvL8LgGOHJQw3suSuYFr+pewHlCFKHCIiabZP579X7w2mJ1wUTTAZoMQhIpINeTDUSDMlDhGRNKhvbPpwYWN99gPJAiUOEZE0WBreGD9xeN/WwlWPBNORn4wgosxR4hARSYMnXnsfgCnHDWktfOGXwfQT34kgosxR4hARSYPHV2wE4JwTD2kt3LQimB5xRgQRZY4Sh4hIGlTtqgVg5MGlQUFixz+zCCLKHCUOEZFMePOZYDr4uGjjyAAlDhGRLqptaPxw4dOzg+kZP8huMFmQ9cRhZoeZ2TNmttrMVprZzLB8gJktMLO14bSDUcJERHLLU6uqADht1MGthe9VBNPR50YQUWZFUeNoAP6Pux8LfBz4ppkdB1wLLHT3UcDCcF5EJOfdW/EuAJ8rPywoaKhrXZhn9zcggsTh7hvdfWn4fRewGjgUOB+YG642F7gg27GJiByI5/93MwBnNT+Ku/h3wfTY8yKKKLMivcdhZiOA8cBLwBB33whBcgEGd7DNVWZWYWYVmzdvzlaoIiKdKi4sCL7M/2EwnTI7umAyKLLEYWZlwIPAd9w96Te5u/scdy939/JBgwZlLkARkSRsq67reGH/w7MXSBZFkjjMrJAgadzt7n8JizeZ2bBw+TCgKorYRERScfuiNwH44slhkljzeDAdeFREEWVeFE9VGXAHsNrd/zth0TxgRvh9BvBItmMTEUnVrc++AcA3JoWJ4v7wZ2z67yKKKPPiERzzVOAyYIWZLQ/LfgDcCNxnZlcC64H8GbxeRPLe0L7FQW/xxrDpanh5tAFlUNYTh7u/AHT0fNrkbMYiItIVb27evW/B0vDB0KEnZD+YLFLPcRGRA/R//7oagBvOC4cVeXRmMP3cHyKKKDuUOEREDtDTa4JneC79+OGwZ1vrggEjI4ooO5Q4REQOwNbdtS3f4wUx+NPFwcwZ/x5RRNmjxCEicgCuefBVAL531tHQ1ASVi4MFp18dYVTZEcVTVSIi3d5Tq4Nmqm+ecRTMD2sZh5+al2NTtaUah4hIil56c2vL9wID/vmbYOYL90YTUJYpcYiIpOjiOf8E4IGvndw6LtXBx0BRnwijyh4lDhGRFGzaWdPyvXx4KfzjlmDmKwsjiij7lDhERFIw9ZfPA/Cz88fAnDOCwuMv7DG1DVDiEBFJ2obte/lgTz0Alx1ZA1UrgwUX3hFhVNmnp6pERJJ0yo1PA3Dj9OPh1hODwovu6hFPUiVSjUNEJAlPrNjY8v2S178dfOlVBmOmRxRRdJQ4REQ6UdvQyNfvXgrA36fXwpvPBgu+/1Z0QUVIiUNEpBPH/fhvAFx8ZAOHPvGloPCK+RDvFWFU0VHiEBHZjxvmraSxyenDHn7x3heDwsnXw0c+Fm1gEVLiEBHpwEPLKrnrxbcpZS8rir8cFB51Fpz2vWgDi5gSh4hIO555vYrv3vsK/djFyuIrg8KhJ8ClD0QbWA5Q4hARaeOxVzfwpTtf5jh7m+XFXw0KD50IX3sh2sByhPpxiIgkuPGJNfz2uTf4WsE8ri28Jygcfxmcf0u0geUQJQ4REaChsYnxP1tAY81u3i6+onXBBbfBuC9EF1gOUuIQkR7v2deruPzOl/lqwaNcV/zn1gXfWw0HHRJdYDlKiUNEeqyqnTV89OcLOTm2kreLZ7cumDADpv06usBynBKHiPQ4H1TXceFtL3L4thd4u/imfRfOWgdlg6IJrJtQ4hCRHmPVhp1cdOtzzOIPPB3/GyR2/P7KM3DohMhi606UOEQkr9XUN/KzR1eyq+Ieft3rN6xs+6v39RdhyJhIYuuulDhEJO9sq67jpseW0+/VO7im8B5mw761i2PPgwt+C0VlEUXYvSlxiEi3V1PfyAOL32T1gju5oukvHBnbyH8AFLau0zToWGIX3QmDj40qzLyhxCEi3Yq7s/SdD/j7M48x+I0HuST+LMXApc0rJIyH0XDsBcSn/BT6H65hMtJIiUNEcpK780bVLpYtfo66VU8wsfp5RsfexYCJ4aftL9juwyZR9qlZMOITYKYfuAzRdRWRyOzcW8fK1avYuvp5YpUvcUT1K4yOvQuAAUeFH+BDI+ttKzua3iddSu/yS6F0IAC6Y5EdShwiklbuzqYPdrF+3Qp2vv0KvmklZTvXcmjtm3wktnmfdQ8CTk4saKc9qap4BA1HnMWQj32WgsM+CrFgpQEZOwPpjBKHiHyIu7N95262bHyLnZveZu+W9TR+UEls13uU7KlkYP1GDvUqCq3xQ9saMDT87KODmwwfxAawpd9Y4iNP5pATP0XRYRMgVtCyfHC6TkrSRolDpJuqq29kx87tVG/fQvWOzdTs3Ebd7m00Vm+B6q3E9m4hXrONoroPKKnfTp+mHfT3HRRZfaf7NqB/+NnvSvuxIT6c7WWjaBp8HGUfOZGhR02gePCR+yQFkjmO5BwlDhGCv7BrG5qoraujrmYP9TXV1NfsoaG2mobaPTTUVNNUv4em2j001e2lsW4PHn6o3wv1e6ChhljDXqxhL7HGGuINNcSbaog31dLLayjyWnp7DSXUtPuXeqp6AYPCT9I6+bFP1ECMLbGD2VE4hL29h9LY5xAKBxxG6ZAjGTj8GPoOHYn1Kulw+0PCj+QfJY4Ubat6j3eWPxPMNDWAO+5NeFMT0ARNTeCNeJPj3oh5MPUmBw+W4R6s403Q5GFZUzDvjjUvc8e8qWVZ8/fmj3lwTGuZb8Raph5OG7FwneblMYLvMW/CCMu8KViP4HvLOrQui7WdT/gUhNvErSnK/zwHzIDi8JPxA6XRXorYZWXsifWhJn4Qdb360tCrHw3FA6BkIPGDBlF00GBK+g+hz4ChHDRwGL2KS8E6DyROB01O0uPlXOIws6nAr4AC4HZ3vzHikPYx4NbjdFOuPWn+QYxKAzHq6EUdhdRaEXVWRJ0VUx8roiFWRH1BMQ2xYpoKimmK98bjvaEw/PQqJVZYTLy4lHhRKYVFJfTqXUZR71KKSkopLulD79J+9OpdBvFenQeThN7hRySbcipxmFkB8BvgLKASeNnM5rn7qmgja7Vyyp+pf+HXOIYTwy0GFgvmzcAKcGJgFiwjXG4WToN5LAYE61gshlsBYGCGWcG+64Ufi7WZb/7eUh4PvxcE7cixAmKxgvAYcYgVBNvECsJ9xbGCGMTimBVgzeWxePi9gFgs2LagIE5BPE48XkCsoJB4QQEF8QLi8V4UFBQQjxcQLygkHo8RLwiPk8RftbkmHn46boARkZxKHMBHgXXu/iaAmd0DnA/kTOIYc8qn4ZRPRx2GiEhkci1xHAq8mzBfCXwscQUzuwq4KpytNbPXshRbrjsY2BJ1EDlC16KVrkUrXYtWx3Rl41xLHO21bfg+M+5zgDkAZlbh7uXZCCzX6Vq00rVopWvRSteilZlVdGX7XBv3qxI4LGF+OLAholhERKQduZY4XgZGmdlIM+sFXALMizgmERFJkFNNVe7eYGbfAv5G8Dju79195X42mZOdyLoFXYtWuhatdC1a6Vq06tK1MHfvfC0REZFQrjVViYhIjlPiEBGRlHTbxGFmU83sdTNbZ2bXRh1PNpnZYWb2jJmtNrOVZjYzLB9gZgvMbG047RGDjppZgZktM7PHwvkeeR0AzKyfmT1gZmvC/z9O7onXw8y+G/7beM3M/mxmxT3pOpjZ782sKrGf2/7O38yuC39LXzezszvbf7dMHAlDk/wLcBzweTM7LtqosqoB+D/ufizwceCb4flfCyx091HAwnC+J5gJrE6Y76nXAYJx3p5099HAWILr0qOuh5kdCnwbKHf34wketLmEnnUd7gKmtilr9/zD345LgDHhNreGv7Ed6paJg4ShSdy9DmgemqRHcPeN7r40/L6L4MfhUIJrMDdcbS5wQTQRZo+ZDQfOAW5PKO5x1wHAzA4CTgfuAHD3OnffTs+8HnGgt5k1Dz22gR50Hdz9eWBbm+KOzv984B53r3X3t4B1BL+xHequiaO9oUkOjSiWSJnZCGA88BIwxN03QpBc6BkvT/sl8H0gcTz3nngdAI4ANgN3hk13t5tZKT3serj7e8DNwHpgI7DD3efTw65DOzo6/5R/T7tr4uh0aJKewMzKgAeB77j7zqjjyTYzOxeocvclUceSI+LABOA2dx8PVJPfzTHtCtvuzwdGErxLqtTMLo02qpyW8u9pd00cPX5oEjMrJEgad7v7X8LiTWY2LFw+DKiKKr4sORWYZmZvEzRXfsrM/kjPuw7NKoFKd38pnH+AIJH0tOtxJvCWu29293rgL8Ap9Lzr0FZH55/y72l3TRw9emgSMzOCduzV7v7fCYvmATPC7zOAR7IdWza5+3XuPtzdRxD8P/C0u19KD7sOzdz9feBdM2se+XQywSsJetr1WA983MxKwn8rkwnuA/a069BWR+c/D7jEzIrMbCQwCli8vx11257jZvZpgvbt5qFJZkccUtaY2SeARcAKWrcO0OwAAASLSURBVNv2f0Bwn+M+4CME/3gucve2N8jykplNAma5+7lmNpCeex3GETwo0At4E/gSwR+IPep6mNlPgIsJnkBcBnwZKKOHXAcz+zMwiWAo+U3A9cDDdHD+ZvbvwBUE1+s77v7EfvffXROHiIhEo7s2VYmISESUOEREJCVKHCIikhIlDhERSYkSh4iIpESJQ3KamU03s+VtPk1m9i8HuL+fmtmZ4ffvmFlJwrLd6Yo7yVhuMLNZKW5zgZn9uINlBxy/md1sZp860O2lZ9HjuNKtmNlVwL8CZ7h7U2frd7KvtwlGUN0Szu9297KuR5n08W8Adrv7zSls8yIwrTnmNssOOH4zOxz4H3efciDbS8+iGod0G2Z2NPBj4LLmpGFmV5vZy2b2atjpCzMbEb6L4n/CdzLMN7Pe4bK7zOyzZvZtgnGMnjGzZxKOMdvMXjGzf5rZkHZi+KiZvRgOIvhicy9tM7vczP5iZk+G7zv4z4RtrjSz/zWzZ8OYbmlnv0eG2y4xs0VmNrqD869NSHQjzewf4fn/rM26H7ouYfmPLHhXxwIL3lMxC8Dd3wEGmtnQ5P+LSE+lxCHdQjg2158IeoevD8umEAyP8FFgHDDRzE4PNxkF/MbdxwDbgQsT9+fuvyYYj+cMdz8jLC4F/unuY4Hnga+0E8oa4PRwEMEfAz9PWDaOoLfyCcDFFrxw6xDgRwTvTTkL+FBCCM0B/s3dJwKzgFvbWedUYGnC/K8IBjQ8CXi/ubCj62Jm5eF1GA98Bihvs/+l4TFE9isedQAiSfoZsNLd70komxJ+loXzZQQ/mOsJBrlbHpYvAUYkcYw64LGEbc5qZ52+wFwzG0UwgmhhwrKF7r4DwMxWAYcTDPnwXMLQDvcDRyfu0IJRjk8B7g+GVgKgqJ1jDyMYNr3ZqbQmxD8Avwi/d3Rd+gCPuPve8LiPttl/FUEtTGS/lDgk54XjUF1IMNLrPouA/3D337VZfwRQm1DUCPRO4lD13nrTr5H2/338DHjG3aeHx3k2YVnbY8Zpf8jqtmLAdncf18l6ewkSV6L2blJ2dF2+28n+i8NjiOyXmqokp1nwboU7gS+GbztM9DfgivAvdszsUDNL5eU8uwj+Ck9FX+C98PvlSay/GPikmfW34G10F7ZdIXyXyltmdhEEox+b2dh29rUaOCph/u8EowJD8MBAs46uywvAeRa8f7uM4M2JiY4GXkOkE0ockuu+RvCmstvaPJJ7cfhWtz8B/zCzFQTvn0glEcwBnki8OZ6E/wT+w8z+TjAy836Fb6P7OcHIxU8RDHO+o51V/xW40sxeAVbS/quQnwfGW2t71kyC982/TEJNpKPr4u4vEwyh/QrBOyoqmmMJ7yEdFZaJ7JcexxXJMDMrc/fdYY3jIYLXADx0gPv6FfCouz/VxVhKCBLRVe6+1MymAxPc/UcHsl/pWVTjEMm8G8xsOUEz0FsE70U4UD8HSjpdq2NzwliWAg+6e/NTWnHgv7qwX+lBVOMQEZGUqMYhIiIpUeIQEZGUKHGIiEhKlDhERCQlShwiIpKS/w9iDU07of7hdQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(solpos['zenith'], pvlib.atmosphere.get_relative_airmass(solpos['zenith'], model='simple'), label='simple')\n",
    "plt.plot(solpos['zenith'], pvlib.atmosphere.get_relative_airmass(solpos['apparent_zenith']), label='default')\n",
    "plt.xlim(0,100)\n",
    "plt.ylim(0,100)\n",
    "plt.xlabel('Zenith angle (deg)')\n",
    "plt.ylabel('Airmass')\n",
    "plt.legend();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
